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Abstract 

We unify and extend a number of approaches related to construct¬ 
ing multivariate Madan-Seneta Variance-Gamma model for option 
pricing. Complementing Grigelionis’ (2007) class, an overarching model 
is derived by subordinating multivariate Brownian motion to a subor- 
dinator from Thorin’s (1977) class of generalised Gamma convolutions. 
Multivariate classes developed by Perez-Abreu and Stelzer (2014), Se- 
meraro (2008) and Guillaume (2013) are submodels. The classes are 
shown to be invariant under Esscher transforms, and quite explicit 
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expressions for canonical measures are obtained, which permit appli¬ 
cations such as option pricing using PIDEs or tree based methodolo¬ 
gies. We illustrate with best-of and worst-of European and American 
options on two assets. 
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1 Introduction 

Madan and Seneta m introduced the univariate “Variance Gamma” {VG) 
process as a model for a financial asset price process with a special view to 
more accurate option pricing on the asset, beyond the standard geometric 
Brownian motion (GBM) model. The VG model has proved to be outstand¬ 
ingly successful in this application, and is in common use by many hnancial 
institutions, as an alternative to the GBM model, despite failing to cater 
for stochastic volatility or a leverage effect, or allowing for temporal depen¬ 
dence of absolute values and/or squares of log returns. Nevertheless, Madan 
and Seneta extended the VG model [H] to a multi-asset version, again with 
a view to important applications in hnance (“rainbow options”), by subor¬ 
dinating a multivariate Brownian motion with a single univariate Gamma 
process (also see [IHl EQl EH [55] )• This construction leaves the marginal 
processes as univariate Variance Gamma processes. But, components are 
dependent by virtue of the common time change. 

Semeraro [Sl| generalises the multi-asset version of Madan and Seneta [H] 
to allow for multivariate subordination. This permits the dependence struc¬ 
ture between asset prices to be modeled in a more flexible way. The economic 
intuition behind multivariate subordination is that each asset may have an 
idiosyncratic risk with its own activity time and a common risk factor, with a 
joint activity time for all assets. In specihc cases it is possible to maintain VG 
processes for each single asset, see jSl] and related applications in Luciano 
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and Semeraro uni SD 112], though this may be sacrihced for more flexible 
dependence modeling, as in Guillaume [25|. Ballotta and Bonhglioni [3| give 
an up-to-date account of modeling based on Levy processes in hnance. 

To summarize, a wide range of multi-asset models based on univariate or 
multivariate Gamma subordination of a Brownian motion has been proposed. 
However, there are still gaps in the literature concerning the characterization 
in general of classes of processes when the class is required to be stable under 
convolution. Further, theoretical results such as formulae for characteristic 
functions. Levy measures and, when possible, transition densities, are needed 
for a comprehensive description of key properties. 

As a prominent link between the real-world and the risk-neutral mea¬ 
sure so as to obtain an equivalent martingale measure, we advocate the Ess- 
cher transform. This approach has been rigorously investigated in the sem¬ 
inal papers of Kallsen and Shiryaev |33] and Eberlein, Papapantoleon and 
Shiryaev im. Esche and Schweizer [18] show that the minimal entropy mar¬ 
tingale measure for a multivariate Levy process preserves the Levy property 
because it corresponds to an Esscher transform. Riischendorf and Wolf EH 
provide explicit necessary conditions for the existence of Esscher parameters 
for multivariate Levy processes and further show that the multivariate Ess¬ 
cher parameter is unique if it exists. Tankov EH provides an introduction to 
the pricing theory in the context of exponential Levy processes. 

The aim of the present paper is to contribute to hlling the above gaps 
by presenting a general class of MAyalued stochastic processes, constructed 
by subordinating multivariate Brownian motion with a subordinator drawn 
from a suitable class of multivariate subordinators. Our intention is to lay 
out a systematic formulation suitable for future development. For the new 
processes, we provide the formulae mentioned in the previous paragraph and 
link the real world and pricing measures by calculating Esscher transforms. 
To illustrate the practical possibilities, we show how the explicit formulae 
can be used to price American and European multi-asset options. 

The most general class of subordinators we consider is Thorin’s [HHl [12] 
class of generalised Gamma convolutions. We call it the GGC class of sub¬ 
ordinators, and the process formed by subordinating Brownian motion in 

with such a process we call a Variance Generalised Gamma convolution 
{y GG) process. Grigelionis [2lj constructed such a VGG-class, which we call 
VGG'^'^ in the present paper. The VGG'^'^ class contains Madan-Seneta’s 
EG as well as the multivariate generalised hyperbolic distributions |1] as 
special cases. Gomplementing Grigelionis’ VGG'^'^ class, we introduce the 
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V class of Levy processes. Our VGG^'^ class includes a variety of pre¬ 
viously derived models such as Semeraro’s a-processes and Guillaume’s 
process [26]. The general VGG = VGG’^’^ U VGG'^''^ class extends the VG 
classes in a number of ways. In particular, the VGG classes allow for inh- 
nite variation and heavy tails. Figure [^depicts the connections between the 
various subordinated classes. 


Var 




VG<^ 


ri 

—► GGGi 
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—► VGG<^'<^ 

VMV% 
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Figure 1: Relations between classes of multivariate Levy-processes. Madan- 
Seneta’s VG occurs as a subclass of Grigelionis VGG‘^'^-class; 

Tf C GGGf Gamma processes and their associated generalised convolution 
processes based on Perez-Abreu and Stelzer Semeraro’s VaV^-class, 
Guillaume’s -class JWf . GF^ = Variance Gamma process based on 

Perez-Abreu and Stelzer E3 ’s multivariate Gamma subordinators, VGG'^’'^- 
class based on based on Perez-Abreu and Stelzer fj^’s GGG^-subordinators; 
V = Variance Matrix Gamma (finitely supported Thorin measures). 
—► points in the direction of generalisation; ■ ■ ■ indicates inclusion in spe¬ 
cial cases. 


The Thorin [581 EH] generalized Gamma convolutions provide a very natural 
class of distributions on which to base our multivariate VG generalizations. 
As we will show, they facilitate construction of a very general class of subor¬ 
dinators and corresponding multivariate Levy processes obtained as subor¬ 
dinated d-dimensional Brownian motions. Our new class complements 1211 , 
and contains a number of currently known versions of multivariate VG pro¬ 
cesses, and extends them signihcantly in a variety of important ways. Al¬ 
though rather technical in appearance, our approach is very much directed 
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toward practical usage of the methodology. Explicit expressions for charac¬ 
teristic functions or Laplace transforms, and Levy measures or densities, are 
derived and exhibited for all our processes. This permits easy programming 
of option pricing routines as we demonstrate by an example, focusing in par¬ 
ticular on the pricing of American style options on a bivariate underlying; a 
thorny problem not often tackled in this context. 

The V GG-class and some of its subclasses are shown to be invariant under 
Esscher transformation, so the risk-neutral distribution constructed as the 
Esscher transformation of a particular member is also in the EGG-class with 
the practical implication that it is possible to calculate prices under the risk- 
neutral measure and estimate model parameters by calibrating to market 
data. Using these concepts, we set up a market model and show how an 
option based on multiple assets may be priced. For illustration we restrict 
ourselves in this respect to a further subclass of the UGG-class which we 
term the UAdTAprocesses. Models in this class have the virtue of allowing a 
quite general dependency structure between the coordinate processes. As an 
example, we price best-of and worst-of European and American put options, 
using a tree-based algorithm. 

The paper is organised as follows. Section contains the theory. In Sub¬ 
section 2.1, we revise the Madan-Seneta VG model and set out two major ex¬ 
tensions: the Variance-Univariate GGC and the Variance-Multivariate GGG 
classes. This necessitates recalling, hrst, some basic facts about Gamma sub- 
ordinators, and then outlining Thorin’s GGC-class. Subsection 2^ gives some 
results on the moments and sample paths of the new processes, including a 
computation of Blumenthal-Getoor indices. In Subsection 2^ we investigate 
relations between our UGG-classes and various multivariate classes of in- 
hnitely divisible distributions, as introduced in Perez-Abreu and Stelzer |15] . 
which we call T^ C GGG^-classes. Subsection 2.4 calculates Esscher trans¬ 
forms, stating the fact that both UGG-classes remain invariant. Subsec¬ 
tion 2.5 introduces the Variance-AdPAsubclass on which we base the option 


pricing model in Section Finally, Subsection 2^ collects further properties 
of our subordinator class, including comparisons with those occurring in the 
literature. 

Section contains applications. Here the market model is introduced, 
risk-neutral valuation is discussed, and in Subsection|3.3|we price some cross¬ 


dependence sensitive options of both European and American types. Some 
illustrations of the kinds of dependencies the models allow are also given 
there. Proofs of Subsections T HH and necessary methodological tools are 
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relegated to Section where polar decomposition of measures and subordi¬ 
nation are briefly covered. 


2 Theory 


2.1 Variance Generalised Gamma Gonvolutions {VGG) 

Preliminaries. is the d-dimensional Euclidean space; elements of are 
column vectors x = (xi,... , 0 ;^)'. Let {x,y) denote the Euclidean product 
and II ■ III; := {x,x) be the Euclidean norm, and set ||a;||| := {x,T,x) for 
x,y G and S G For A C let = 74 \{ 0 }. 1a denotes the 

indicator function. The Dirac measure with total mass at x G is 6x- 
/ : M —)■ M denotes the identity function. 

X = {Xi,..., Xd)' = {X{t))t>o is a d-dimensional Levy process if X 
has independent and stationary increments, X( 0 ) = 0 and the sample paths 
t ^ X(t) G are cadlag functions. 

The law of a Levy process X is determined by its characteristic function 
via = exp{f'0x(d)} with Levy exponent, for t > 0, 9 E 

i’x(e) = n,f(di). (2.1) 

Here 7 x G M'^, Ex G is a symmetric nonnegative dehnite matrix, Hx is 
a nonnegative Borel measure on satisfying 



A 1 nx(dx) < oo , 


( 2 . 2 ) 


and II ■ II is a given norm on We write X ~ L'^( 7 x, Ex, Hx) when¬ 
ever X is a d-dimensional Levy process with canonical triplet (yx, Ex, Hx); 
BM'^Gi •= 0) refers to Brownian motion with drift 7 and covari¬ 

ance matrix E. 

Paths of X are of (locally) finite variation (FV^) whenever Ex = 0 and 

I ||x|| nx(dx) < 00 . (2.3) 

In this case, we write X Fl/'^(Dx,nx) with Dx denoting the drift of 
X: Dx ■= 'yx — /o<|p||<i^ Llx(dx) G M'^. A d-dimensional Levy process 


6 


T with nondecreasing components is called a d-dimensional subordinator^ 
possibly with drift Dt, written T ~ A general Levy process 

X ~ L'^( 7 x, Sx, Ifx) is a subordinator with drift Dx if and only if X ~ 
FV‘^{Dx, Ifx) with Dx G [ 0 , ooY and fix being concentrated on [ 0 , oo)^. 

We write X = Y and X ~ Q whenever C{X) = C{Y) and C{X) = Q, 
respectively, where C{X) denotes the law of a random variable or stochastic 
process X. There is a correspondence between inhnitely divisible distribu¬ 
tions and Levy processes X: for all f > 0 the law of X(t), P{X(t) G dx), 
is inhnitely divisible. Vice versa, any inhnitely divisible Borel probability 
measure Q on determines uniquely the distribution of a Levy process via 
X(l) Q. This connection is used throughout the paper. For instance, 
we write T ~ Qs to indicate that T is a subordinator with T(l) ~ Q. 
See [21 m [121 [13 EHl [52] for basic properties of Levy processes and their ap¬ 
plications in hnance. 


Subordination. In [ 6 ] various kinds of subordination are introduced (see 
Subsection 4.2 for details). In the present paper, we will make use of two 
extreme cases: univariate and (strictly) multivariate subordination. Let 
X = (Xi,..., Xd)' be a d-dimensional Levy process. X serves as the subor¬ 
dinate. 

Given a univariate subordinator T, independent of X, dehne a d-dimen- 
sional Levy process, denoted X o T, by setting 


(X o T)(f) := (Xi(T(f)),..., X,(T(f)))', t>0. (2.4) 

In the sequel, we denote the law of X o T by C{X) o £(T). We refer to this 
type as univariate subordination (cf. Section 6 in [52]). 

Suppose X has independent components Xi,..., X^. Let T = (Ti,..., T^)' 
be a d-dimensional subordinator, independent of X, and dehne a d-dimensional 
Levy process by setting 

X O, T:=(Xiori,...,X,or,)'. (2.5) 

The law of X T is denoted by C{X) jC,{T). 

Remark 2.1. When dealing with strictly multivariate subordination, we 
have to restrict the class of admissible subordinates X to Levy processes 
with independent components. This is necessary if we are to stay in the class 
of Levy processes. For instance, let B ~ BMi{0, 1) be a univariate standard 
BM. Then X = (B, B)' is a Levy process, but 1 1 —>■ {B{t), B{2t))' is not. □ 
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Gamma subordinator. Denote by r(a,/3) a Gamma distribution with 
parameters a,P > 0, i.e., a Borel probability measure on M, absolutely con¬ 
tinuous with respect to Lebesgue measure, dehned by 

r(a,/5)(da;) = l(o,oo)(a:) xeR. (2.6) 

We write G ~ r5(a,/9) for a Gamma process G = {G{t))t>o with parameters 
a, P > 0, that is, G is a univariate subordinator having marginal distributions 
G{t) ~ r(at,/9), t > 0. Further, for A > —f3, f > 0, recall that 

= {(/3/(/3 +A)}“* = exp|-f / (l - (2.7) 

V 0 

(the hrst formula is well known, the second identity is Frullani’s integral, 
see p.l 6 & p.73 in [ 8 ]). In particular, a Gamma process has zero drift, and 
its Levy measure admits a Lebesgue density nG(dr) = l(o,oo)(’^) ae ^^drjr. 

For a = (3 we have F^[G(1)] = 1. G is then also called a standard Gamma 
process, briefly G ~ F 5 (a) := F 5 (a,a). 

Madan-Seneta VG^ Process. Madan and Seneta |11] (for extensive in¬ 
vestigations and reviews cf. Ha EDI EH EH EH sal E5]) suggest subordinating 
Brownian motion with a Gamma process. For the parameters of this model 
we assume pi G M'^, b > 0 and S G with S being symmetric and 

nonnegative dehnite. 

Let B ~ S) be a d-dimensional Brownian motion and G ~ F 5 ( 6 ) 

be independent of B. A Levy process D is a d-dimensional Variance Gamma 
(GG‘^)-process with parameters b, pi, S whenever Y = B o G, which we write 
as 

Y ~ VG'^ib, /i, S) := BM'^iia, S) o Fs(6). (2.8) 

(a) Note that a GGAprocess has zero drift and is of finite variation. 

(b) The Laplace transformation of Y takes on an explicit form, straight¬ 
forwardly derived from conditioning: for f>0, AgM'^ with |||A||| —(/i. A) < b, 

Eexp{-(A,F(t))} = | 6 /( 6 +(;,,A)--||A|||)} . (2.9) 

(c) If S is invertible, for f > 0 the distribution of Y[t) is absolutely continuous 



with respect to the Lebesgue measure with the following density: 

2M/H“ exp{{S-V.!/)} f ll!/ll|-. 

mt){y) 7r'^/2 (detS)V2 r(6t) \ 2 b+\\^\\l_J 

xAV-,i|/2(v^(26+l|Al||-.) Nil-.), !/€«■'. (2.10) 

Further, still with det E 7 ^ 0, the canonical Levy measure of Y is absolutely 
continuous with respect to Lebesgue measure and satisfies: 

dny _ 6 2<"-'‘>/" exp{(S-‘fty)} [• 26+11^111-, 

dy ^ 7r‘^/2 (detS)V2 | ||2/|||-i j 

X AV2(v^(26+||fi|||-.) I|!/|l|-.), S/eR" (2.11) 


Here is the modified Bessel function of the second kind (see |23], their 
Equations (3.471)-9 and (8.469)-3; also |T5], their Appendix). It is conve¬ 
nient to revise the following facts about the modihed Bessel function Ki, of 
the second kind: for 5 ,7 > 0 , i/ G M 

POO 

2{6/'yy^‘^ exp { - (5/r) - yr} dr . ( 2 . 12 ) 

Jo 

For large values of x (see |23], their Equation (8.451)-6), 

Ky{x) ^ Kxi 2 {x) = \pn e~^/\f2x , x —)■ cx). (2-13) 

We use a variant of Ky^ defined by 

Ky{r) -.= r'^Ky{r), r, z/> 0. (2.14) 


Remark 2.2. We get from (2.7) that the shifted Gamma process G{t) = 
at -|- G{t), t > 0, with G ~ F 5 (q;, /3) has Laplace transform for A > 0 


— log Ee = Xat + at log = Xat + t 

P 


qOjOo) 


log(l + {X/y))a50{dy ). 


For independent processes Gi,...,G„ with Gk = T s{oik, fik)) 1 < k < n, 
introduce the associated shifted Gamma processes by Gkifi) = Okt + Gkif), 
t >0, l<k<n. It follows from the independence that for A > 0 


— logEe = — ^logEe = tXa+t / log(l-|-(A/j/)) 7)j(dj/), 

k=l J{0,oo) 
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where a = J2k=i ~ Sfc=i discrete measure associated 

with the increasing function Un{y) = Ylk=i^ki[i3k,oo){.y), 2 / > 0. Using the 
representation in the last display and taking suitable limits in distribution, we 
arrive at the class of generalised Gamma convolutions. In the multivariate 
setting in Subsections 2.5 2.6, we shall employ hnitely supported Thorin 
measures, such as 7^ above, in the context of the A^F'^-class. □ 


Generalised Gamma Convolution Subordinator (GGC). The class of 
Gamma distributions is not closed under convolutions. To extend the Madan- 
Seneta UG'^-class it is convenient to use the subordinators corresponding to 
Thorin’s class [SHI EH] of generalised Gamma convolutions {GGC). This is 
the smallest class of distributions that contains all Gamma distributions, but 
is closed under convolution and weak convergence (see uniEiEiiianiEaiES]; 
for multivariate extensions see |5iniE9]). The class of GGG-distributions is 
a subclass of the self-decomposable distributions and, thus, the distributions 
are inhnitely divisible. 

A d-dimensional Thorin measure T is a Borel measure on [0, oo)^ with 



(l+log ||xi|) A (l/||x||) r(dx) 


< oo. 


(2.15) 


(Throughout x = x'^ — x~ denotes the decomposition of an extended real 
number x G M into positive and negative parts.) 

A subordinator T is a GGG^-subordinator with parameters a and T, in 
brief T ~ GGGg{a,T), when T is a d-dimensional Thorin measure, a G 
[0, oo)'^ and, for alH > 0, A G [0, ooY, 


logUexp{-(A,T(f))} = t{a,\)-\-t 




|x|p -I- (A, x) 


\x\ 


■os ■! II 111','’’ !• T(Ax) . 

(2.16) 

The distribution of a Thorin subordinator is determined by a and T. Any 
Thorin measure T admits a polar representation T = a K, relative to 
:= {x G [0,cx))‘^ : ||x|| = 1} (see Lemma 4.1 below). This allows us 
to specify the corresponding Levy measure. For T GGGl{a,T) with 
T = a (8)p /C we have T ~ ^^(a, IIt) with 


dlfTo — 
k{s,r) = 


'Si Jo 


dr 

6rsk{s,r) — a(ds), (2.17) 


' (OjCxd) 


e /C(s, dr), r>0,sG§+, (2-18) 
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see isms], their Theorem F and Proposition 4.3, respectively. 

Variance—Univariate GGC (VGG^’^). As a hrst extension of the VG^- 
model, we review Grigelionis • El class. Grigelionis used univariate sub¬ 
ordination o and subordinated Brownian motion with a univariate GGC- 
subordinator. For the parameters of his model we take /i G a > 0 and 
S G with S being symmetric and nonnegative dehnite. Further, let 

T be a univariate Thorin measure. Let B ~ S) be a d-dimensional 

Brownian motion, T GGGg{a,T), independent of B. 

Given such B and T, we call a Levy process of the form Y = B o T a d- 
dimensional Variance Univariate Generalised Gamma Convolution {VGG'^’^)- 
process with parameters a, T. We write this as 

U ~ /i, S, T) := S) o T). (2.19) 


The next theorem gives the characteristic function and Levy density. Part (a) 


(our Equation (2.21)) is proved in Subsection |4.3[ Part (b) occurs in [21j (see 
his Proposition 3). (In (2.20), log : C\(—cxd, 0] —)■ C denotes the principal 
branch of the logarithm.) 


Theorem 2.1. Let Y ~ UGG^’i(a, p, S, T). 

(a) For all 9 G t>0, we have Eexp{i {9,Y{t)) = exp{t-^y (6*)} with 


^l)Y{0) = ia {fr,9) - (a/2)||0||| 



r-i(/i,g) + |||g||| 

T 


r(dr). (2.20) 


(b) Assume det S 0 and T 7 ^ 0. Then fly is absolutely continuous with 
respect to d-dimensional Lebesgue measure on , where, for y G Mf, 


'^(y) = 21 ^-')/^ (detE)->/^ ||j,||^^7 exp{(E-V.9>} x (2.21) 

X [ (2r+||,*|||.0"'‘/fi/4v'(2r+|H|_.)ll!/ll|-.)r(dr). 

J{0,oo) ^ 

Besides this, we have fly = {J^ ^rs 5 'y('S, r) dr/r) Q!d,E(ds) in Euclidean 

polar coordinates with Lebesgue surface measure ad^E on (for d = 1 we 
interpret ad,E os the counting measure), for r > 0, s G 

^y(s,r) = ^-d /2 (detS)-i /2 ||^||-d^ exp{r (S”V, s)} x (2.22) 

X [ .Ad/2((2r+||/i|||_i)^/2 ^||s||e-i] r(dr). 

J(0,oo) ^ ^ 
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Remark 2.3. In both classes, VG‘^ and VGC^’^, we subordinate a Brow¬ 
nian motion with a single univariate subordinator. Thus the components 
of these processes must jump simultaneously. To allow the components to 
jump independently of each other we must use multivariate subordination 
of Brownian motion. This motivates our next step, the introduction of our 
VGG^’^-c\ass. □ 


Variance—Multivariate GGC {VGG^^^). Next we give another modihca- 
tion of the VG'^-model which is constructed by multivariate subordination o^. 
For the parameters of this model we assume a d-dimensional Thorin measure 


T, /i G 


e [0, 00)“^ and S G 


'bdxd 


with S = diag(Sii,..., Y^dd) having 


nonnegative entries. (We impose on Brownian motion the requirement to 
have independent components, so as to stay in the class of Levy processes. 


see Remark 2.1 
Let B 


be a Brownian motion. Let T ~ GGGg{a,T) be 
independent of B. Given such B and T, we call a Levy process of the form 

D 

Y = B Od T a d-dimensional Variance Multivariate Generalised Gamma 
Convolution (VGG^’^)-process with parameters We write this as 


Y ~ VGG^’^a, /r, E, T) := BM^fi, S) GGGl{a, T). 


(2.23) 


To state formulae for the characteristics of this process, dehne the outer 
o-products oi y = (?/i,..., yd)\ z = {zi ,..., ZdY G and S G as 

yoz := {yiZi,y 2 Z 2 ,...,ydZdyeW^, (2.24) 

E 02 ; := diag( 2;i,..., ^d)E G . 

We can decompose [0, cx))^ = U 0 ^jc{i d} semi-cones, where 

Gj := I Xjtj : > 0 for all j G j| , 0 7 ^ J C {1,..., d} , (2.25) 

and Cj are the unit coordinate vectors. Analogously, we can decompose = 
U 07 jc{i,...,d} Vj into Vj := ■ Xj ^ 0 for all j G J} for 0 ^ J C 

{1,..., d}. Let be the cardinality of J. 

We need a family of reference measures. With £ denoting the univariate 
Lebesgue measure dehne £j := <S)k=i as the product measure with the 
following factors: 


£j,k := lAk)£ + lii,...,d}\j{k)6o , 1 < < d. (2.26) 
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Observe £j(Mf — Vj) = 0. Finally, to provide an analog of (2.22) in Euclidean 
polar coordinates, let aj^E denote the Lebesgue surface measure on OVj 
for J 7 ^ 0. (In the discrete case we interpret aj^E as counting measure). 

The next theorem gives the characteristic function of Y and an expression 


for its Levy measure. It is proved in Subsection 4.3 
Theorem 2.2. Let Y ~ /r, S, T). 

(a) For all 6 G t>0, we have Eexp{i {6,Y{t))} = exp{t^jJY{0)} with 


'iPy{0) = i(hoa,^)-2ll^llsoa- 


[0,oo)^ 


log-[ 731^5 -^- ' (da;) 


a; 


(2.27) 

(b) Always, fly = Y.$^jc{i,...,d}^J’ nj(Mf - Vj) = 0 for ^ ^ J C 
{!,..., d}. If T{Cj) = 0 then flj = 0. Otherwise, if T{Cj) > 0 and 
det E > 0 then flj is absolutely continuous with respect to £j and, for y G Vj, 


dn. 


fy) = 2^^ #0/2^ #J /2 exp{(E V;!/)} x 

f Tidx) 


(2.28) 


r n 


2 ||x||^ + (/iOa:,E V) I 


#a/4 


X 


K#j/2 ( I (^2||a;||2 + (/ioa;,E ^ 


j&J 


1 / 2 ^ 


whereas in flj = /sdnVj(/o°° dj(’S, r) dr/r)aj^E{ds), /or r > 0, s G S|nVj, 


gj{s,r) = 2(2-#0/2 ^-#J/2 exp{r (E'V, s>} x 
f T{dx) ( n,, _ ,'I-#a/2 


(2.29) 


n /2 1/2 


'Cj UjeJ^/j iej 




X 


1 y ^2 

K#j/2{r X |(2i|x|p + (/iOx,E-V>) ) • 


3&J 


2.2 Moments and Sample Paths 


In Subsection |2.6| , we see that both EGG-classes support pure jump processes 
with inhnite variation and inhnite moments. In Propositions |2.1|j2.2 


we 


provide conditions on the Thorin measure that can be used to check local 
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integrability of fl'p and Ily as well as existence of moments (for a proof see 
Snbsection 4.4). Set || ■ ||i := | ■ | and || ■ ||rf := || ■ ||. 

Proposition 2.1. Let t > 0, 0<q<l, p > 0 and T ~ GGCg{a, T). Then: 
(^) /o<||z||<i ll^ll'^ nr(dz) < cx) /\\xp < oo. 

(b) .E[||T(t)f] < CX) /Q^||^||<^r(dx)/||crf < cx. 

Proposition 2.2. Let kG{l,d} andY r~^VGG‘^’’^{a, p,T,T). Then we have: 
(a) Let 0<q'<2. If 

[ r(dcr)/||x||^/^ < X 




then 


'0<||j/||<l 


Ill/ll'' ny(d|/) < X. 


(2.30) 

(2.31) 


If, in addition, det S > 0 then (2.30) and (2.31) are equivalent, 
(b) Let p,t>0. If 

f[o i]d < °° when p = 0 

T(da;)/||a;||^ < x when p ^ 0 


(2.32) 


then 


E[r(f)r]<x 


(2.33) 


If, in addition, either Y\k=iLk 7 ^ 0, or p = 0 and det S > 0, then (2.32) 
and (2.33) are equivalent. 


Remark 2.4. Snbordinating Brownian motion with independent snbordina- 
tors from the generalised inverse Ganssian (GIG)-snbordinators one obtains 
the class of generalised hyperbolic Levy processes (0 El DEI for detailed 
acconnts). Halgren |27] identihed the nnivariate GIG distribntions as gener¬ 
alised Gamma convolntions, with the implication that the associated class of 
hyperbolic Levy processes 0 forms a snb class of the RGG'^’^-class (see [TO] 
and [21|, his Example 1). 

Owing to Proposition [2^ we may restrict onr analysis to the snbordinator 
class. For {a, /3, 7 ) G M x (0, x)^ U (0, x)^ x {0} U (—x, 0) x {0} x (0, x) 
the GIG{a, ft, 7 )-distribntion has the following probability density on (0, x) 


GIG{a, ft, 7 )(dx) 




«,/3,7 


X 


"“-1 da ;, 


cr > 0. 
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Here Ca,p,'y is a normalising constant. For 7 = 0 the identity GIG{a, (3, 0) = 
r(a,/9) holds, and exponential moments are hnite for A < /3, and 

this extends to other parameters a ,/?,7 as long as /? > 0. For (3 = 0 and 
a < 0 observe GIG{a, Oj'y) =invF(— 0 , 7 ) is the inverse Gamma distribution 
with hnite p-moments of order p < \a\ only. 

We determine the Blumenthal-Getoor index [5] as an indicator for activity 
of the associated G/G-subordinator, for potential applications see p m\. 
In (2.7), note x'^nG(da;) = a is hnite for all 0 < g < 1, 

forcing the Blumenthal-Getoor index of GIG{a, (3, 0) = F(a, (3) to degenerate 
to zero. For the remaining parameters, where 7 > 0, we compute the Laplace 
transform of GIG{a, (3,^) for A > 0 as 

exp{-A„,^,^(A)} := / e"^"=G/G(a,/3,7)(da;) = — 

Jo ^a,/3-|-A,7 


= 27“/2G«,^, 


iLH(2V7(/5 + A)) 
' (/5 + A )“/2 


x/^G « 7 .( 2 a-i )/4 exp{- 2 V 7 (/j + A)} 

V7rO„,/3,T,7 ^ ;^y2a+l)/4 ’ 


A —)■ 00 , 


as follows from (2.12)-(2.13). Observe AQ,_y 3 ^g,(A) ~ 2 a/^ as A —)■ 00 and 
thus cxd)) rsj 2y^7/x as a; 4 , 0 by a Tauberian 

theorem (see |H], P-75) with associated Blumenthal-Getoor index 1/2. □ 

Remark 2.5. Other possible extensions of the RG^-class comprise a range of 
possible sample path behaviour. In |15] the univariate GGMR-processes have 
been identihed to be subordinated Brownian motions, and it is also known 
that the associated subordinator is a GGG^-subordinator (see [31], their 
Example 8.2). As perceived in |T3|, Blumenthal-Getoor indices of GGMY- 
processes exhaust the whole of the interval (0,2). In particular, for any 
given q G (0,1) there are GGG-subordinators with Blumenthal-Getoor index 
g, the latter by Part (a) of Proposition 2.2| (see H for multivariate GGMY- 
models.) □ 


Remark 2.6. With (|2.15|) being straightforwardly verihed. Too is a Thorin 

dx 


measure on (0, cxd), where 


770(dx) := l(o,i/e)(a;) 


dx 


xlog'^(l/x) 


Y l(e,oo) ( 2 ^) 


log (x) 
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As T^{dx)/x‘^ = Too{<^x)/xP = oo for 0 < g < 1 and p > 0, re¬ 

spectively, we get from Proposition 2.1, that any associated GGC^{a^Too)- 


subordinator has inhnite p-moments, with Blumenthal-Getoor index equalling 
1. By Proposition 2.2, also any 0, 0, S, 7^)-process with det S > 0 

has inhnite p-moments, and its Blumenthal-Getoor index equals 2. □ 


2.3 Levy Process Classes via Polar Decomposition 


Based on the polar decomposition of the Levy measure (see Lemma 4.1), 
Perez-Abreu and Stelzer |19] construct classes of self-decomposable distribu¬ 


tions on cones generating classes of Levy processes of hnite variation [FV^), 
surpassing subordinators, and including versions of Gamma and GGG pro¬ 
cesses. 

Multivariate Gamma Process. We reproduce the model in [49] (their 
Section 3). Recall := {x G : ||x| 


= 1}, and let /?:§“—)■ (0, cxd) be a 


Borel function, and a a hnite Borel measure on such that 


1 -I- /3(s) ^ 


< 


oo 


(2.34) 


We refer to a d-dimensional Levy process W as a pAprocess with parameters 
a and /?, written as W ~ P^(a,/3), whenever X ~ Fl/'^(0,nx) with 


By = 



dr 




Srs — a(ds) 


(2.35) 


Assuming a,/3 satisfying (2.34) it is shown in [49] that the RHS in (2.35) 


dehnes a Levy measure (see their Proposition 3.3). Much in the spirit of our 
the P^-class carries processes with inhnite moments, amongst 


2.2 


Subsection 

other things (see H. their Examples 3.14, 3.15 and 3.16). 

Multivariate Gamma Subordinator. Recall := [0, ooY fl and let 

/3 : —)■ (0, cxd) be a Borel function, and a a hnite Borel measure on 

such that ( 2.34[ ) is satished, but with replaced by 

We refer to a d-dimensional subordinator T as a pAsubordinator with 
parameters a and /?, written as T ~ Tg{a,(3), whenever, for all A G [0, oo)*^. 


-logEe-AGG = f 


/3(s) -I- (A, s) 
log— 


(2.36) 
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As follows from the Frullani identity in (2.7), the RHS in (2.36) matches 


(3{s) + (A, s) 

log—= 


dr 


'Si Jo 


a(ds), 


and, thus, T 


Tg{a,(3) holds if and only if T ~ with 11^ as 


in (2.35), but with S“ and If^ replaced by §1 and If-r, respectively. 


Plainly, a P^-subordinator is a P^-process with nondecreasing compo¬ 
nents: S"^nP^ = P^. In the univariate case, we have P5(a<5i, (5) = P 5 (a, id(l)). 
The connection with our GGC'-class is P|(a,/9) = GGCg{0, a <5/3(.)). 

The associated lAPf-class is thus a subclass of our RGG'^’Adass. Comple¬ 
menting this, it is possible to contrive the RP^-class as a class of Levy pro¬ 
cesses associated with the matrix-gamma-normal class of [19] (see their Sub¬ 
section 5.3.2). 

Multivariate Gamma Convolution Process. In [lU] (see their Dehni- 
tion 4.4 and their Subsection 4.4) it is also shown how to characterise the 
associated class of multivariate generalised Gamma convolutions associated 
with cones. It is possible to introduce a GGG|-class of FPApi-ocesses, ex¬ 
tending the polar decomposition from [0, oo)!^ to , as we illustrated in con¬ 
text of the inclusion P^ C P^. In particular, we have fl GGGf = GGGg. 

We conclude this subsection by establishing further inclusions in Figure 
postponing further investigations of this kind to Subsection |2.6 
Variance Gamma (VG) revisited. Assume Y ~ VG^{b, such that 
Y ~ VGG'^’i(0,/i,S,r) with T = b6b in ([2T^. If d = 1 and S ^ 0 


then (2.22) degenerates to, for r > 0, s G = {±1}, 


gY{s,r) = 5exp |r(s/i — (26S-I-/i^)^'^^|s|)/s|y|s|, 


and in (2.35) we have 

^l{b (d_i + di), ±1 ^ ((2&S + p2)i/2 ^ , 


Y 


comparable to the representation of univariate VG^-processes as a difference 
of independent Gamma subordinators in 


For d > 1 and invertible S we get from (2.13) and (2.22) that, for s G 
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as r —)■ cx), 


X 


S-i 


(2vr)('i-i)/2(detS)V2 


(d+l)/2 


Plls-i 


X exp |r( (S V,s> - (2&+||/i|||-i)^/^ l|s||s-i)} 


The asymptotic equivalence in the last display does not match (2.35), erasing 
the possibility of Y being a T^-process. 

Within his l/GG'^’^-class Grigelionis (see |23], his Proposition 3) shows 
that any VGC^'^iQ, p, S, T)-process is self-decomposable, provided either d = 
1 or d > 2 and /i = 0. Assuming an invertible S and imposing further 
moment conditions upon T, which are duly satisfied for finitely supported 
Thorin measures such as T = bSb, he shows that his result is sharp: if /i 7 ^ 0 
and d> 2 then a GGG'^’^( 0 ,/i, E, T)-pro cess cannot be self-decomposable. 
Consequently, a /i, S)-process with d > 2,/x 7 ^ 0 and invertible S 

cannot be self-decomposable, let alone be an element of the P^ C GGCf- 
classes. Plainly, iPi /2 is completely monotone. By differentiating (2.12) under 


the integral sign we find that iP^'^ 2 ( 0 +) < 0 - particular, if d > 2 , then 

Kd /2 is no longer completely monotone, and a CCAprocesses with d > 2 , 
/i = 0 and invertible S cannot be a GGC^-process. 

It would be interesting to provide a detailed study regarding the VGG'^’^ 
{k G {1, d})-classes and the P^ C GCC'^-classes, but this is beyond the scope 
of our present paper. 


2.4 Exponential Moments and Esscher Transformation 

We use the notation, for one hence all f > 0, 

= {A G < cx)} = {A G : f < 00 }. 

d|hll>i 

(2.37) 

Vy is a convex subset of M'^, containing the origin (see p. 165), similarly 
we introduce Vy. 

Further, we need to introduce 

C>A := { 2 : e [0, 00 )^ : ||a:||^ > (A, x)} , A G , (2.38) 
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and a transformation 


<Sx{x) 



X, 


X eRi. 


(2.39) 


We provide conditions on the Thorin measnre ensnring finiteness of exponen¬ 
tial moments for the associated GGC/VGG-mode\ (see Snbsection 4.5 for a 
proof). 


Proposition 2.3. Let t > 0, k E A G T ~ GGGi(a,T) and 

Y r~.VGG'^'^a,fi,E,r). Then: 

(a) {0} U ([0, cxD)f\(9A) is a convex and compact subset o /and S\ is a 
continuous function from Ox into [0, cx))f. 

(b) A G Vt -v^ simultaneously, T([0, cxD)f\(9A) = 0 and 



log ||5a(x)|| r(da:) 




T{dx) < oo. 


(2.40) 


(c) For A G Vt the image measure of the restriction T\Ox{-) '■= T{Ox H ■) 
under the mapping Sx, denoted by Tx '■= (T|C>a) o is a well-defined 
Thorin measure on [0,cxD)f. 

(d) Without restrictions on (a,p., S, T): A G Vy qx,k £ Vt, where 


f + il|A|||, ifk = l, 
\ A o p-1-4S(A o A), if k = d 


Assnme that Y ~ L'^^'yy, Sy, fly) is a Levy process with respect to an 
nnderlying stochastic basis (ff, V, {Vt}, P). 

The Esscher transform on V with respect to Y is given by 


Y 


(^Qx,t 

dP 


exp (A,y(t)) 
Epexp (A,Y(t)} ’ 


For t > 0 and A G Vy in (2.37) it is well-known that 

[ 0 , 1 ]. 


t > 0 , A eVy. (2.42) 

Vt ^ [0,1] 

dehnes a probability measnre, eqnivalent to P : Vt —>■ [0,1]. Besides this, 
{y(s) : 0 < s < t} remains a Levy process nnder the new measnre Qa*- 
Next we show that both l/GG-classes are invariant nnder Esscher trans¬ 
formations (for a proof see Snbsection |4.5 ) We provide more specihc ex¬ 
amples thronghont the remaining part of the paper, more specihcally see 
Theorem Snbsection 2^ and Section 
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Theorem 2.3. Lett>0, kE{l,d}. IfY r^VGG‘^' 
then q := qx^k £ and {T(s) : 0 < s <t}\Q\^ ~ 
with qx^k £ and Tq as in Proposition 


2.3. 


^(a, yU, S, T) with X G Vy 
l^G'G"^’^(a,/i + SA,S,7;) 


2.5 FXr^-Class 

In this subsection we restrict ourselves to finitely supported Thorin measures 
and consider a corresponding subclass of VGG^’'^. 

A^F'^-Subordinator. Let n G N = {1,2,...}, 6* = ( 6 i,..., bn)' G (0, oo)"-, 
M = {mki)i<k<d,i<i<n e having columns Mi ,..., Mn e [ 0 , cx))f and let 
Gi ~ r 5 ( 6 i),... ,Gn ~ r 5 ( 6 „) be independent standard Gamma processes. 

We call a d-dimensional subordinator T an -subordinator with pa¬ 
rameters n, 6 *, M, briefly T rsj Air|( 6 *,M), provided 

n 

T ^ M(Gi,..., G„)' = GiMi . (2.43) 

1=1 


Next, we show that AlF'^-subordinators are GGG'^-subordinators, but having 
zero drift a = 0 and hnitely supported Thorin measure: 


Lemma 2.1. Let T ~ M) with n = dimfe*. Then T ~ 

GGGg{0,7T), where, simultaneously, 

^‘'(o,nT) = 

n p 

ffr = exp{-6;r}dr/r, 

J{0,ca) 

n 

(2.44) 

Tt = '^bl5^,^Ml/\\M^\\^ , 

l = l 

(2.45) 

n 

Vt = fjlAGM'": {Mi,X)<bi} 

i=i 

(2.46) 

and, for t > 0, A G Vy, 


n 

- logi7e<"'^W> = tY,h log{(6; - {Mi,X))/bi] . 
1=1 

(2.47) 


Proof. By (2.16), (2.47) follows straightforwardly from (2.45). 
of Proposition (4.1) is applicable to (2.43) giving Ilr = X]r=i ^ 


Part (a) 
G,Mn with 
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a similar superposition and intersection in place for 7 t and Vt, respectively. 
It remains to verify the formulae in (2.44)-(2.46), but with n = l, M = Mi, 
b = bi and G = Gi ~ ^sib). 

With n ;= 6 /(Qgo) dr/r and A G with (A,M) < b calculate 


(l_e(A,x>)n(dx) = b / e-'"’dr/r 


(2.48) 


On the other hand, i^exp (A, G(f)M) = i? exp{(A, M) G(f)} in which we 
can substitute the characteristic exponents of G using (2.7). By (2.48), the 
resulting expressions match those in (2.44)- (2.45). □ 

For the remaining part we review some properties of the Gamma distribution. 


Lemma 2.2. Let c, a, ai, ..., a„, Pi,. .., Pn > 0. Let Z ~ r(a, P). Let 
Zi,... ,Zn be independent with Z^ ~ L{ak, Pk) for 1 < k < n. Then we have 
cZ ~ V{a,P/c) as well as the equivalence ‘(i)^(ii)’, where: (i) YfJk=i^k ~ 
T{a,b) for some a,b > 0; (ii) Pi = ■ ■ ■ = Pn- If (i) or (ii) is satisfied then 
b = Pi and a = Yf'k=i “fc- 


Proof. Note aSb = Ylk=i holds if and only if both b = Pi = ... = Pn 

and a = 'Yl'k=i^ki and the equivalence ’(i)-v^(ii)’ follows from the Thorin 
representation of the Gamma distribution. □ 


We introduced the Tg{a,P) in (2.36), and its connection with the GGG- 
class was Tg{a,P) = GGGg{0,a Z)p Si 3 (.)) with the identity r; 5 (Q; 5 i,/ 3 ) = 
r 5 (Q!,/ 3 (l)) in the univariate case (d = 1 ). 

As we clarify next, a given AdF^-subordinator does not need to have 
Gamma marginals nor does it need to be F^-subordinator. 


Lemma 2.3. Let T = (Ti,..., Td)' ~ AiTg{b^, M) with n = dim 6 *. 


(a) Then (i) (ii), where (i) Tk ~ Ts{j)k, Qk) for some pk, qk e (0, cx)); 

(ii) there exists l<lo<n with mk^i^ > 0 such that bimk^i^ = bi^mk^i for all 
l<l<n with mk,i 7 ^ 0 . 

(b) (P) ^ (ii’), where (P) T r-, Tg{a,P) for some a,P; 

(ii’) for all l<k,l<n, the following implication holds 


\\Mi\\Mk = \\Mk\\Mi => \\Mi\\bk = \\Mk\\bi. (2.49) 


In addition, if one of (i) or (ii) holds then we have qk = bi^/mk^iQ and pk = 
"d/so, if one of (P) or (ii’) is satisfied then a = JP(=i biSMi/HMiH 
and’p{Mi/\\Mi\\)=bi/\\Mi\\ {l<l<n). 
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Proof, (a) follows from Lemma 2.2 as we can decompose the kth. component 
of T in (2.43) into a sum of n univariate Gamma subordinators. 

(b) ‘(ii’)^(i’)’: For S ~ r|(a,/3) and x G [0, with Euclidean norm 
ll^lll: = {x,x) = 1, introduce a univariate subordinator by 


S^{t) := l{„,:„>o}(A^(s)) (a;,A^(s)) , f>0. 


(2.50) 


0<s<t 


To determine the Levy measure fl^ of S^, let A C (0, cx)) be a Borel set and 
note 

rL^(y4) = E[jf{0 <t<l: AS{t) G {ax : a > 0} and (x, AS'(t)) G A}] 

= n 5 ({r G [0, oo)^ : r G {ax : a > 0} and {x, r) G A}) . 

As I| a3;:a>0|(r5) = l{„x:a>o}(s) = l{x/||x||}(s) for r > 0, s G we get from 
(2.35) that the RHS in the last display matches 

r r°° Aj, 

n^(R) = / / lA(r’(x,s)) l{^/|p||}(s) 

Jsf Jo ^ 

poo 

= a({x/||x||}) / l^(r/||x||) dr/r . 


Substituting r = r'||x|| on the RHS of the last display, it follows from (2.7) 
that either = 0 or 5® ~ r 5 (Q!({x/||x||}), ||x||/?(x/||x||)). 

Thus prepared, let T ~ A^rg(fo,M) with T = (2.43) for 

independent standard Gamma subordinators Gi,... ,Gn. Let T* be dehned 

)y T. In particular, 

Gi{Mk,Mi)/\\Mk\\E, l<k<n. 


as in (2.50), but with S replaced by T. In particular, 

j^Mk/\\Mk\\E P 


E 

In addition, suppose T ~ Tg{a, f3). Then must either be degener¬ 

ate or a univariate Gamma subordinator. Gonsequently, by Lemma |2.2 


we 


must have bi\\Mk\\% = bk {Mk, Mi) for 1 < / < n with ||Mi||Mfc = ||Mfc||Mi. 
The latter is equivalent to (ii’), completing the proof of ‘(ii’)=»(i’)’. The 
proof of ‘(i’)=^(ii’)’ is analogous. □ 

Variance—A4r {V With parameters fo = (6i,...,6„)' G (0, cxd)” 
and M G as set for an WlTAsubordinator, in addition, take /i = 
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(yUi,..., Hd)' e and a diagonal matrix E = diag(Sii,..., Yidd) with non¬ 
negative entries. 

Whenever Y = BodT, with B, T being independent and B S) 

being Brownian motion, while T ~ M), we call Y a Variance MY 

( VMr^)-process^ written in the following as 

Y ~ VMV\K, M, /i, E) := BM\^i, E) M). (2.51) 

For a generic case, where det E 7 ^ 0, we give formulae for the canonical Levy 
measure fly. To each column Mi we associate both a dimension l<di<dhy 

di := #{1 < k < d : ruki > 0} , 1 < ^ < n, 

and a cr-£nite Borel measure := <S)k=i^’ik ^ product measure 

with the following factors 


tii, := 1(0,oo){mki)iYl{o}{mki)So, l<k<d,l<^<n. (2.52) 

For 1 < / < n, we set 

I3i := 2bi+ ^ rUkipllYlkk = 2biY{ino Mi,E~^n) , 

a, := 6,/3,(2-53) 


The next theorem gives formulae for the Levy measure and Laplace exponent 
of Y, which has hnite variation (recall (2.3)) and is invariant in form under 
Esscher tranformations. 


Theorem 2.4. Assume Y Vpi,YYj with n = dimb^. Then: 

(a) We have Y ~ /r,E,T) with T = YA=lhlS^,^Ml/\\Ml\\'^■ 

(h) Always, Y ~ Fl/'^(0,ny). When, in addition, det(E) > 0, then for all 
Borel sets A C 


ny(Ll) = 

n 

1=1 


^di/2 


Pi 




dk! P^kk^kl)^ 


di/4 


exp 


^ ^ikVk/Y^kk \ ti{dy). 


23 





(c) We have 


Vy = {X e : (/ioMi, A) + -\\X\\y^m^ <bu 1 < / < n} , 
and, for t>0 and A G Vy, 

n . 

- logEe^^’^W) = t^6jog|(6i-(/ioMi, A)--||A|||^A^J/6i|, (2.54) 

i=i 


and 


{F(s) : 0<s<t}|Q;,^ 




Here nx = /i+SA, and Mx G [0, oo)'^^"' has the following columns Mf ..., M^: 


Mf = 


bi (/ioMi,A) 2ll'^llsoM; 


M, 


1 < I < n. 


(2.55) 


Proof, (a) follows from Lemma 2.1 and (2.51). AlL'^-subordinators T have 
zero drift. Also, /o<||a;||<i ||a;||^/^n 7 ’(da;) is finite for AdLAsnbordinators T. 
Thus, Y ~ Fl/(0,ny) by Part (c) of Lemma 4.2, (A stronger result follows 
from Part (a) in Proposition 2.2 ) In view of Part (a), the remaining parts (b)- 
(c) follow from Theorems 2.2 2.3 as well as Proposition 2.3 □ 


Remark 2.7. If R ~ I/A1P‘^(6*, M, p, E) with n = dim6* then 


y = J2^i, (2.56) 

1=1 


as an implication of (2.54). Here Yi, 
VG<^{bi,^ioMi,EoMi) 


, Yn are independent processes with 
< n. It is, thus, possible to con- 


Yi ~ VG‘^{bi,fvo Ml, E o M;) for 1 < / s n. 
struct RAdPAprocesses by superimposing independent Madan-Seneta VG^ 
processes. Wang [6T] comes to similar conclusions, and constructs multivari¬ 
ate Levy processes with RG^-components by superimposing suitable RG^ 
processes, just as in the right hand-side of (2.56). □ 


Transition Densities. For a subclass of AdPAdass it is possible to obtain 
formulae for transition densities, as we illustrate next. 
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Let 6* = (6i,..., 6rf+i)' e (0, cx))'^+^ and M G [0, such that, 

simultaneously, 

M = {mki)i<k<d,i<i<d+i = (diag(mii,... ,mdd),Mrf+i) and 

d 

Wrukk mk(^d+i) Y 0- (2-57) 

k=l 

With t > 0 dehne 

d 

c; -.= c;(b„M) -.= {6j«Vr(*j+i)}n''h/(r(«'>‘)™h). (2-88) 

k=l 

d 

(3* := /3*{K,M) := -bd+i + y^Yk^k,d+i/mkk- (2.59) 

k=l 


The proof of the next result follows from a similar analysis as in Section 48.3.1 
in [36]. (We are unable to provide substantial simplihcation of the integral 
in ( 2.60[ ) which occurs by integrating the joint density of d + l independent 
Gamma random variables. However, using the results in [IH], it is possible 
to expand the integral in terms of Lauricella functions.) 


Lemma 2.4. Let t > 0 and T ~ JLiTg{b^, M) with M satisfying (2.57). 
Then T(t) admits a Lebesgue density fx: for r = (ti, ..., r^)' G 


d 

fmir) = C; l(o,oo)‘^(^)exp { -^bkTk/mkk]y- 

k=l 
d 

gtfed+i-i Y[{Tk-mk{d+i)sy’"^~^ ds . (2.60) 

k=l 




Next, we state such formulae for the associated model. Let /i = 

(/ii,..., fid)' and E = diag(Eii,..., Tdd) be the parameters of the underlying 
Brownian motion. To ensure the existence of a Lebesgue transition density 
for the HAlT'^-process, we need to make the additional assumption that E 
is invertible, i.e. all E^^. > 0. 

With the help of ( |2.58| ) and (2.59) dehne 


Ok ■ ^ / (y^Oik^d+lf'^kk) ) 

Ck ■ 26/j T TlikkLk/^kk ; 


Ok := mk(d+i) ((5fc/nrfcfc) + (/ifc/(2Efcfc)), 
Cfc . 'Y Ck/(Tikkriikk) ■ 


25 










Further, for t > 0, we set 


Ct := C';(6 „M) ( 2 VdetS)-i/ 2 j^ 


m 


tbk-\ 

k{d+l) 


k=l 


k=l 


Cd+1 ■— 


d 

2&(i+l “1“ ^ ^ 'l^k{d+l)f^k/^kk ) 


k=l 


Dk := 2 n c 


-d (d-2ba+it)/4 


d-\-l 


d+1 


n ('>i‘Vr(M)) 


X 


k=l 


X 


n 

k=l 


Jl-26fct)/4 (3+26fet)/4 (l+2fefet)/4 -1/2 \ 


'‘kk 


kk 


m 


Theorem 2.5. AssumeY ~'FA^r'^( 6 *, M,/i, E) wth det S 7 ^ 0 and (2.57). 
Then the law ofY{t) admits a Lebesgue density fyit) fort>0: 


fy({)(.y) = Qgxp (y^A‘tyt/Sn;}x 


(2.61) 


k=l 




d „1 


n 

k=l ' 


ex' 


akuyl afcs) (1 -m) 


\tbu — l 


u J u^^k +2 


-du ds 


A exp { y^/ifcr/fc/Sfcfcj 


X 


(2.62) 


k=l 


K\2ba+^t-d\/2{\JCd+l Yjk=l {^kkmki^d+l))) ^ iF|26fet-l|/2 {CklVk “ ^fc|) dZk 


/or?/ = (?/i,...,?/rf)' e 


n 

k=l 




Proof. Since we have Y = B o^T for independent B BM'^{/i,'E) and 
T ~ A^r'^(&*M), (2.61) follows from Lemma 2.4 by conditioning Y{t) = 
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B{T{ t)) o n values of T{t) (see Part (b) of Lemma 4.2). In view of Re 

V 


mark 


2.7 


we may write Y = Yi + Y 2 for independent processes Yi and Y 2 . 
Here the d-dimensional process Yi has independent components with the kth 
component being a /ifcmfcfc, Sfcmfcfc)-process (l<fc<d). Further, Y 2 

is a VC^i^bd+i, p. o Md+i,T, o Md+i)-pTocess. The formula in (2.62) follows 
from ( 2 . 10 ) by convolution. □ 


2.6 Subclasses of GGC'-Subordinators 

In this subsection we review subordinator classes as they occur in the lit¬ 
erature and relate them to our formulations. The GGC^, T'^, and 
subordinator classes were introduced in Subsectionsand 2.5[ respectively, 
and one of the conn ectio ns was Tg{a,(d) = GGGg{0,a < 5 / 3 (.)). As dehned 
in (2.43), in Lemma 2.1 the AlTAdass was identihed to be the subclass of 


GGCAsubordinators with drift a = 0, having hnitely supported Thorin mea¬ 
sures T. As we clarihed in Part (b) of Lemma 2.3 a given AdPAsubordinator 


does not need to be a PAsubordinator. Two other classes, such as the ones 
introduced by Semeraro and Guillaume [26] , are related to them as shown 
in Figure]^ (Compare Figurej^with Figure[^) In the univariate case, where 
d = 1, note that aF^ = F;^ = F 5 . 


^ MY^s 

aTi —► ^F^ : GGG^ 

^ rf ^ 

Figure 2: An arrow points in the direction of generalisations of different 
subordinator classes, as described in the text. ■ ■ ■ indicates inclusion in special 
cases. 


Perez-Abreu and Stelzer’s F^-Subordinator 

nators based on 


The Fg-class of subordi- 


was dehned in (2.36). In the univariate case, we have 
observed F 5 (a< 5 i,/ 3 ) = Fs'(a,/d(l)). The connection with our GGG-class was 
F|(a,/3)=GGG|(0,a®p5^(.)). 

Let T ~ FKa, /9), A G M'^, q = q\ & T>t as in ( |2.41 ). We get from Part (a) 
of Proposition 2.3 that, simultaneously, afff) < ((?,•)}) = 0 and (2.34) 
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holds with /3 replaced by /3x{-) := (3{-) — (g, •). In Part (c) of Proposition 2.3 
observe that = a^pSjS^. Consequently, the associated CP^-class 

of subordinated Brownian motions is closed under the Esscher transform in 


the interpretation of Theorem 2.3 


Semeraro’s CK-Subordinator. Semeraro [m] introduced another approach 
to multivariate Gamma subordinators (also see [1011111112]). The parameters 
of this model are as follows: let a, 6 G (0, cxd), = (ai,..., ad)' G (0, oo)*^ 
such that, simultaneously, b > aau for all 1 < k < d. Let Si,, Sd+i be 
independent such that 


~ Psl —-a, — ), l<k<d, 
'Oik ak' 


Sd+i ~ rs(a,6). 


We refer to T as an a-subordinator, in brief T ~ ar|(a,6 ,a*), provided 


T = (Ti,...,Trf)' with 

Tk := Sk + akSd+i. (2.63) 

Any aPAsubordinator T admits standard Gamma marginal distributions: 
Tk ~ ^s{b/ak). As a result, the associated GaP^-processes, called aVG in 
have GG^-marginal distributions. 


We give an alternative representation of T in (2.63). Introduce para¬ 


meters 6* = {bi,... ,bd+iy G (0, and independent standard Gamma 

subordinators Gi,... ,Grf+i, Gk ~ ^s{bk) for l<k<d+l, by setting 


bk := -a, 1 < k < d, 

Oik 


bd+l ■ ® ) 


and, with Si,..., Sd+i as above, 

b 


Gk ■— 


b — aak 


Sk , 1 < k < d, 


Gd+i := 


d+l 


For T in (2.63) we conclude that T ~ Adr|(6*, M(a, 6, a*)), where in our 
notation 

M{a,b,aS) := (^ia.g{b—aai,...,b—aad),aa^/b G [0, . (2.64) 

We show that the V process is not closed under Esscher transform 
by considering the following bivariate example. In Part (c) of Theorem 2.4 
we have A := (1,0)' G Vy for p = (0, 0)', S := diag(l, 1) and 

T~ar|(l,2,(l,l)')=V(rj((l,l,l)',(^ 1 / 2 ))' 








but also, recalling (2.55), 


— (M; 


- M. 


(i,oy ^(i,oy 


2/3 0 2/3 \ 

0 1/2 2/3 ) ■ 


T) 

Under the Escher transform, the hrst component of T reads Ti = (2/3)Gi + 
(2/3)G3 for independent Gi^G^, ~ r 5 (l) by (2.43). However, as implied by 
Lemma 2.2, (2/3)Gi + {2/‘i)G^ ~ r5(2,3/2), and Ti cannot be a standard 
Gamma process. The associated Uar|-class of subordinated Brownian mo¬ 
tions is thus not closed under the Esscher transformation in the interpretation 
of Theorem 12.31 

Guillaume’s Subordinator. Guillaume [20] extends Semeraro’s aT'^-class 
as follows: let a* = (ui,..., a^)', a* = (ai, • • •, a^)',= (/^i, ■ ■ ■,/3d)' G 
(0, oo)'^. Cl, C 2 > 0. Let S'!,..., S'd+i be independent such that 

S'fc ~ r5(afc,/3fc), l<k<d, S'd+i ~ r5(ci,C2). 

We refer to T as a -subordinator, in brief T ~ ^rg(Q;*, a*,/3*, ci, C 2 ), 

provided T = (Ti,... ,Td)' with := Sk + UfeSd+i- 

With Si,..., Sd+i as above, introduce independent standard Gamma sub- 
ordinators Gi ~ rs(ai),..., Gd ~ r 5 (ad), Gd+i -^Ts{ci) by setting 


_/3fe o 

v_j^ ' 

dk 


1 < k < d, 


Gd+i ■— —Sd+i 

Cl 


We conclude that T ~ M(a*, a*,/3*, ci, C2)), where in our notation, 

5* = (oi ,... ,ad, Cl)' e (0, cx))'^+^ and 

M(a*,a*,/3*,Ci,C2) := (^diag (ai//3i,..., ad//3d) , (ci/c2)a*j G [0, _ 

(2.65) 

Further, observe that 

{QTs{a^,a,,,/3^,Ci,C2) : a^,a^,(3^ e (0, 00 )'^, Ci, C 2 > 0} 

= {Air5(6*,diag(a;*),2/*)) : x*, y* G (0, cx))'^, 6* G (0, cx))'^+^} . 

By Part (c) of Theorem 2.4, the U^T'^-class of subordinated Brownian mo¬ 


tions is, thus, closed under the Esscher transformation in Theorem 2.3 


Unlike the aT'^-subordinator, a given ^T'^-subordinator does not need to 
have Gamma marginals, and we clarihed this in Part (a) of Lemma 2.3 By 


Part (b) of Lemma |2.3[ contemplating (2.65), a ^P'^-subordinator and, thus, 
any aP'^-subordinator is also a Pg-subordinator, concluding settlement of our 
diagram in Figure]^ also recalling the chain of inclusions aP'^ C ^P^ C WIP^. 
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3 Applications 


We are primarily concerned with demonstrating how our 1/Wir'^-subclass 
can be applied, in particular, to price multi-asset options. The V 
subclass, as we showed, contains other popular models, such as the multi¬ 
variate VG m, the Semeraro aVG la, and Guillaume’s extension ra. 

In Subsection ^ a market model using the GWlT'^^-process is introduced, 
and we give explicit formulae for the expected value of the fc-dimensional 
log-price process and its covariance matrix, and for the expected value of 
the price process itself. This allows us to tabulate values of these quantities 
for a specihc parameter set which we will use to illustrate the results. The 
corresponding densities are calculated using the formula for the characteris¬ 
tic function given in (2.27) of Theorem 2.2 and displayed in Figure]^ The 
parameters required to make the Esscher transform an equivalent martin¬ 
gale measure linking the real world and risk neutral dynamics are derived in 
Proposition |3.2 of Subsection 3.2 As an example, pricing of four kinds of 
two-asset options, specihcally, European and American best-of and worst-of 


put options, can then be operationalised as we demonstrate in Subsection 3.3 


The exact form of the Levy measure as given in Theorem 2.4| (b) is an essen¬ 
tial ingredient here. 


3.1 A WAir^-Market Model 

We employ the GAdreprocess to model the log-prices of risky assets of a 
hnancial market. Potentially latent risk factors are described by a process 
Y ~ M,/i, E)-process with respect to a given stochastic basis 

(r2,W, {Wt},P). The risk factors drive a fc-dimensional price process S with 
Siit) = S'i(O) e^**'*^ for f > 0 and i = l,...,/c, with /c-dimensional log-price 
process R = [m — q + k)I -|- AY = {R{t) : t > 0} given by 


R{t) = {m — q + k) t + AY(t) = {m — q + K)t + X{t), t >0 , (3.1) 


where m G M*’ is the expected total return rate of the assets, q E is the 
dividend yield of the assets, A G with rows A* G M'’* satisfying A^' G Vy, 
i = determines the factor loading of the corresponding log-return 

process, and k G is an adjustment vector given by Ki = — logPe^**-^^ = 
— logPe^^* such that ESi{t) = S'j(O) f > 0, i = 1,..., k. Recall 

/ : [0, oo) —)• [0, cx)) denotes the identity function. Proposition 3.1 gives 
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formulae for the moments of R{t) and the explicit form of the adjustment 
vector K. 


Remark 3.1. The dependence structure of the risk factor process Y is lim¬ 
ited, as S has to be a diagonal matrix in order that we remain in the class of 
Levy processes. The matrix A maps those risk factors to specihc asset prices 
and generates a richer and perhaps more realistic dependence structure, for 
similar arguments and setup see [IQISTIEII. Accordingly, AY and R are not 
necessarily RAdT^^-processes, but are of course Levy processes. □ 


Proposition 3.1. Let t >0 and R as in (3.1) with n = dim 6 *. Then: 

(a) ER{t) = {m — q + K + A ^1^=1 T ^ ^i) t > 0. 

(b) Coy{R{t)) = A fcr=i U{ToMi){fioMiy + ToMi 


A't, t>0. 

(c) Ki = YJi=ih log{(&i - (/ioMi,A*') - l^i] ) * = 


Proof. This follows by differentiating the Laplace transform (see (2.54) in 
Theorem 2.4). □ 


p 

ERiil) 

ER^il) 

Var(Ri(l))^ 

Var(R2(l))^ 

Cor(Ri(l),R2(l)) 

0.30 

0.0917 

0.0782 

0.1296 

0.2104 

0.3651 

0.00 

0.0921 

0.0780 

0.1260 

0.2114 

0.0329 

-0.30 

0.0919 

0.0785 

0.1276 

0.2092 

-0.3076 


Table 1: Expected value, volatility and correlation of R{1) for A = 
(1, p; p, 1)°'®, p E {—0.30,0,0.30}, Y ~ M, p, E) with parame¬ 

ters n = 3, d = k = 2, m = (0.1,0.1), q = (0,0), 6 * = (5,5,10)', 
M = (0.5, 0,0.5; 0,0.5,0.5), pi = (-0.14,-0.25), S = diag(0.0144, 0.04). 


We investigate the distribution of R for parameters: d = k = 2, m = 
(0.1, 0.1), q = (0, 0), 6* = (5, 5,10)', M = (0.5, 0, 0.5; 0, 0.5, 0.5), p = (-0.14, 
-0.25), S = diag(0.0144, 0.04) and A = ( 1 , p; p, l)^-^ with p e {-0.3, 0, 0.3}. 
Table states the expected value, volatility (square root of variance), and 
correlation of R(l), for p G {—0.3,0, 0.3}. These numbers facilitate a bet¬ 
ter understanding of the potentially abstract model parameters and serve as 
a basis for comparison when the Esscher transform is discussed in Subsec¬ 
tion 3.2 The expected values for both coordinates are below m = (0.1, 0.1) 


and are robust when varying p. The expected value of the hrst coordinate be¬ 
comes maximal for p = 0 whereas for the second coordinate the relationship 
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is inverted. This effect is determined by the term A ^11=1 h ^ Proposi¬ 


tion 3.1 (a). A similar behavior can be observed for the volatilities, however, 


here the roles of the coordinates are exchanged. Most notably, the correla¬ 
tion differs considerably from the dependence parameter p. The main driver 


of this difference is the first component A ^ ^ 


A' 


m 


Proposition 3.1 (b). Depending on the sign of the entries of Ap this term 
increases or decreases the correlation. For p G {—0.30,0, 0.30}, Ap has nega¬ 
tive entries in both coordinates, consequently increasing the correlation above 
p. This effect weakens when decreasing the dependence parameter p. 

Figure illustrates the density of R for t G (0.01, 0.25} when varying 
p G {—0.30,0,0.30}. The densities are obtained numerically from the char¬ 


acteristic function given in (2.27) of Theorem 2.2 using fast Fourier inversion 


For t = 0.01, the superposed processes ApoT dominate AYA^'^B o^T, where 
T ~ AlFg(6*,M) and B is d-dimensional standard Brownian motion. For 
p = 0, most of the probability mass is located near the x- and p-axes. For 
p = 0.30, additionally mass appears around two straight lines in the first and 
third quadrants (positive dependence). For p = —0.30, additionally mass ap¬ 
pears around two straight lines in the second and fourth quadrants (negative 
dependence). For t = 0.25, the density is close to normal with nearly ellipti¬ 
cal level lines. Note, though, that for p = 0 the density is not symmetric but 
skewed towards the left and lower values. 

Remark 3.2. A desirable property of a parametrisation of a multivariate dis¬ 
tribution is to be able to distinguish between parameters describing marginal 
distributions, and parameters describing the dependence. For the V 
however, this is in general not possible. Each parameter appears in at least 
one marginal distribution. This is a consequence of the fact that the family 
of Gamma distributions is not stable under convolution, except for singular 

for 

□ 


cases; see Lemma 2.2, These are the cases analysed by 0. See also 
correlating Levy process and related applications. 


3.2 Risk-Neutral Valuation via Esscher Transform 

Option pricing requires a risk-neutral measure as the basis for risk-neutral 
valuation. In the general Levy process setting, such a measure is not guar¬ 
anteed to exist and further, if it exists it is in general not unique. But in 
Part (c) of Theorem 2.4 we showed that the RAtFAdass is invariant under an 
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Figure 3: Density level lines of R(t) = {m — q + K)t + AY(t) for t G 
{0.01, 0.25}, A = (1, p; p, 1)0-5, p ^ {-0.30, 0,0.30}, Y ~ p, S) 

with parameters the same as for Table [I[ 


Esscher transformation, and here we follow common practice by adopting the 
Esscher transformation for identifying a risk-neutral measure, see [I51E21EI]. 

For the processes R,X,Y in (3.1) and h G Vr = Vx = Day fhe Esscher 
transform is given by (see (2.42[)) 


dQ^ 


h,t 


Ph,R{t)) 


Ph,X{t)) 




dP Ep[eAMt))] Ep[eX^d))] Ep[eA'h,Y{t)^ ’ 
such that, with h G Dp = Dx = D^Y) 


f > 0, 


dg« dg 

dP 


X 

h^t 


dg 


Y 

A'h^t 


dP 


dP 


for f > 0 


By Part (c) of Theorem 2.4, as Dp = Dx = Day, we observe that 
Vn ={h e K‘ : {^oMi,A'h} + tp'/illL,,, < k, l<l<n 


(3.2) 


(3.3) 


Also, by replacing A with A'h in Theorem 2.4, it follows from (3.3) that 
{Y{s):0<s<t}\Ql, ~ EMF"(&,^M^/i^E'^), heDp,t>0, 
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1 < / < n. 


with 6!^ = 6*, = /i + S A'h, = S, and 


Mt = 


bi - {fioMi,A'h) - l\\A'h\\l^^^ 


Ml 


Next, we summarise risk-neutral pricing under the Esscher transform: 

Proposition 3.2. Assume h* E such that h*,Ci + h* E Vji = Vx = 
T^ay, 1 < i < k. Then, for the market with price process Sq = and 
Si = Si{0) with S'j(O) E M"*", I < i < k, the Esscher transform is an 
equivalent martingale measure with respect to the numeraire Sq: -r ~ P 

and e^^^Si/ Sq are rp-martingales, for 1 <i < k and T > 0 if and only if 


mi - r = AAY{i){ti) + AAY{i){h*) - AAY{i){ti +h'^), forl<i<k, (3.4) 


where Ax is the cumulant-generating function of an W^-valued random vari¬ 
able X, i.e. Ax{u) = u E {v eM.^^ : < oq}. 


Proof. Let h E Vay such that h Ci E Vay, for l<i< k. Then ■= Qhx 
is well-dehned and P( 3 ^|e^**S'i(f)/S'o(f)| <oo, for l<i<k and 0<f<T. Note 
that e^^^Si/SQ is the exponential of a Levy process, under both P and Q/ft, 
and thus for 1 < i < k and 0 < f < T it is the case that 


Eq, [e^^^S,{T)/SQ{T)\Xi] 


e'^'^Si{t) 

Soit) 

e^^^Siit) 

Soit) 



^*(l)/5.(0)' 

5o(l)/5o(0). 


T-t 



Ep [e<^N(i)>] 


S„{t) y Bp |e<'‘."-(i))] j 

6 '^ Si{t) (jy,.j^,.,_._„.j^AAY{n+h)-AAY{h)){T-t) 

SQ{t) 


T-t 


Recall Ki = — logPe^^* 
e^^^S,/SQ is a Qh-martingale, i 


— logPe^'^'P^*^^^^ = —Apiy(ej) to see that 
1,..., k, if and only if h satishes (3.4). □ 


Remark 3.3. The parameter h* is called the Esscher parameter. For general 
exponential Levy market models. Theorems 4.4-4.5 in [33] ( also see |3I], their 
Theorem 2.6) state that h* is unique, provided the driving Levy process does 


34 













not degenerate under P in the sense of Definition 24.16 of [52] • An application 


of this result yields that our market model (3.1) admits a unique h*, provided 
rank(y4) > k, rank(M) > d and det S > 0. 


□ 


Next we set the interest rate to r = 0.05 and keep the remaining model pa- 

The resulting Esscher parameter, the adjusted 


rameters as in Subsection 3.1 


risk-neutral parameters and some basic statistics are provided in TableThe 
Erst row indicates the three different scenarios, i.e. p G {—0.30,0,0.30}. In 
the second row the Esscher parameter h* is seen to have negative components 
that are increasing in p. The sign of the components of h* is as expected, 
since the model under P corresponds to a “bullish” market with expected 
return rates of irti = m 2 = 0.1 exceeding the risk free rate r = 0.05, and h* 
has to counterbalance this effect. The third row gives the transformed pa¬ 
rameter p!^ which tends to be lower than the original parameter under P and 
is increasing in p as well. The matrix distributing the Gamma subordinators 
to the coordinates is displayed in the fourth row. The elements are all 
greater than those of M and the more negative the dependence parameter 
p becomes the stronger is this effect. The resulting characteristics of the 
distribution are displayed in Rows 5 to 8. These numbers can be compared 
to the numbers under P in Table The expected values of R(l) under the 
Esscher martingale measure are lower than under P. The volatilities increase 
across the board by nearly 1%. For the correlation the same can be observed; 
an increase of about 1% is found when comparing the Esscher numbers to 
the original numbers under P. Summarising, volatilities and correlations in- 



Table 2: Esscher parameter and resulting basic statistics for A = 
(1, p; p, 1)°'®, p G {—0.30,0,0.30}, r = 0.05, Y ~ VM, p^Y) with 
parameters the same as for Table 
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crease when we change from P to Q^*. Thus under the pricing measure Q^* 
risk in the form of volatilities requires a higher risk premium than would be 
anticipated under P, e.g., when pricing a call or put option. Further, diversi- 
hcation effects are less pronounced under the pricing measure, e.g., requiring 
a higher premium for basket options. 


3.3 Pricing Best-of and Worst of Put-Options 


The hnancial market model presented above can capture a wide range of 
dependencies between different asset prices. As an illustration we price some 
cross-dependence sensitive options of both European and American styles. 
European options can be conveniently priced by Fourier methods na. Thus, 


we can draw on the results provided in Theorem |2.2| to compute European 
option prices. Pricing American options can be carried out by hnite differ¬ 
ence methods, discretising the respective pricing partial integro-differential 
equations, or by using tree-based methods. See |2H] for a recent survey on nu¬ 
merical methods in exponential Levy process models. Both methods require 


formulae for the Levy measure that we provided in Theorem 2.4 


As an example we consider best/worst-of put options with respective early 
exercise values 


xboij,i(«) = ( -R' - V = ( R - A 


(3.5) 


2=1 


2 = 1 


for 0 < f < T, where T is the maturity date and K G M"*" the exercise price. 

The risk-neutral parameters are: n = 3, d = /c = 2, 6* = (5,5,10)', 
M = (0.5, 0, 0.5; 0, 0.5,0.5), S = diag(0.0144,0.04), /i = (-0.14, -0.25), m = 
(0.1, 0.1), q = (0, 0) and A = (1, p; p, 1)°'® with p G {—0.3, 0, 0.3}. Note that 


we have set here r = 0.1 in contrast to Subsection 3^, resulting in A* = 0 and 
Qh* _ p allows us to interpret the option price dependencies on the 

parameter p without confounding this with effects of the Esscher transform 
on the option premium. To compute American option prices we use the tree 
approach as outlined in [32l |3l], based on [16]. The European option prices 
are obtained as a byproduct of this procedure. 

The recombining multinomial tree calculation we use has probability 
weights derived from the Levy measure, as provided in Theorem |2.4[ The 
option parameters are set to T = 0.25 and K G (90, 95,100,105,110}. The 
tree models the bivariate process Y = (1^,12)^ directly, with an exponential 
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p 

K 

Best-of put price 
European American 

Worst-of put price 
European American 

0.3 

90 

0.04 

0.05 

0.75 

0.81 

0.3 

95 

0.18 

0.24 

1.76 

1.90 

0.3 

100 

0.71 

1.06 

3.74 

4.03 

0.3 

105 

2.17 

5.00 

7.00 

7.49 

0.3 

110 

4.98 

10.00 

11.32 

11.98 

0 

90 

0.01 

0.02 

0.76 

0.82 

0 

95 

0.09 

0.13 

1.83 

1.98 

0 

100 

0.44 

0.77 

3.96 

4.27 

0 

105 

1.63 

5.00 

7.48 

7.96 

0 

110 

4.27 

10.00 

12.01 

12.62 

-0.3 

90 

0.00 

0.01 

0.77 

0.83 

-0.3 

95 

0.03 

0.06 

1.85 

2.01 

-0.3 

100 

0.24 

0.53 

4.14 

4.45 

-0.3 

105 

1.19 

5.00 

7.94 

8.42 

-0.3 

110 

3.66 

10.00 

12.63 

13.20 


Table 3: Best-of and worst-of put option prices for T = 0.25, K G 
{90,95,100,105,110}, A = (l,p;p,l)°•^ p G {-0.30,0,0.30}, r = 0.10, 
Y ~ M,/i, S) with parameters the same as for ra6/e[^ 
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transform to obtain the price process. At each node of the tree the process 
branches on a regnlar rectangnlar 127 x 127 grid. The minimnm step sizes 
are 4.92 x 10“^ and 8.37 x 10“^ for Yi and I 2 respectively. Prices are then 
obtained to an accnracy of three signihcant digits. The time increment is 
1.25 X 10“^. Rnn times are rednced by trnncating propagation of the tree in 
its spatial dimensions after one time increment. Allowing the tree to grow 
fnrther does not affect the resnlts. 

The resnlts are presented in Table As expected, pnt options prices 
are increasing in the exercise price K. Also, the worst-of pnt option prices 
exceed the corresponding best-of pnt option prices, which is consistent with 
no-arbitrage. For ont-of-the-money options, the early exercise preminm is 
higher for the worst-of pnt compared to the best-of pnt. The early exercise 
preminm for at-the-money options is approximately similar in both cases. 
For in-the-money options, the early exercise preminm is higher for the best- 
of pnt compared to the worst-of pnt. The dependence parameter p affects 
the option prices as expected. The payoff of the best-of pnt increases the 
contingency that both price processes fall jointly, thus the option premium is 
increasing in p. The payoff of the worst-of put increases if at least one price 
process falls, thus the option premium is decreasing in p. 

4 Proofs 

4.1 Polar Decomposition of Measures 

For /i, 1 / being a-hnite measures, p® v denotes the corresponding unique a- 
hnite product measure. The trace held of the d-dimensional Borel held 
in A e is denoted by SJ4, and = {x G : ||x|| = 1} is the unit 

sphere for a given norm || ■ || on We say that a Borel measure p is locally 
finite relative to B ^ provided p{C) < 00 for all compact subsets 

C C B. Let /C : X -BJqoo) [0; 00 ] b® ^ locally finite Borel transition 

kernel relative to (0, 00): simultaneously, s h-)■ }C{s,B) is Borel measurable; 
B I—)■ /C(s, B) is a Borel measure, locally hnite relative to (0, 00). There exists 
a unique measure a ® )C : ® i3((0, cxd)) —)■ [0, cxd], locally hnite relative 

to X ( 0 ,cxd) satisfying {a®lC){A x B) = J^IC{s, B) a{ds) for A G 
B G i3((0,oo)) (for example see Exercise 3.24, Chapter III of |29]). Dehne 
a®plC B(R.f) —)■ [0, 00] as the image oi a® 1C under x (0, cxd) 9 (s,r) 
rs G By construction, a ®p JC is a locally hnite Borel measure relative 


38 


to satisfing J^dJ{x) {a K:){dx) = f{rs) K:{s,dr)a{ds) for 

nonnegative Borel functions /. 

Next we provide a polar decomposition on as a disintegration of 0 ^ 
for Borel measures satisfying additional integrability conditions. The result 
is directly applicable to Levy and Thorin measures, as in Lemma [4T] we may 


choose w{r) = r^Al and w{r) = (l+log“ r)A(l/r) in view of (2.2) and ( |2.15[ ), 
respectively. We omit the proof. It is possible to adapt the arguments in [5] , 
Lemma 2.1, and [50], Proposition 4.2, respectively. 


Lemma 4.1. Assume 0 < /jgdta(||x||)/i(da:) < oo for a Borel measure fi on 
and a continuous function w : (0, cxd) —)■ (0, oo). Then we have: 

(a) fi is locally finite relative to with /r(Mf) G (0,oo]. 

(b) There exists a pair {a, (3) such that, simultaneously, 

(i) a is a finite Borel measure on 

(a) 1C : X i3((0, oo)) is a Borel kernel, locally finite relative to (0, oo); 

(Hi) 0 < f w(r) /C(s, dr) < oo for all s G 
(iv) fi = a ®p 1C. 

(c) If {a',1C') is another pair, simultaneously satisfying (i)-(iv), then there 

exists a Borel function c : —)■ (0, oo) such that a(ds) = c(s)Q;'(ds) and 

c{s)lC{s,dr) = lC'{s,dr). 


4.2 Subordination and Decomposition 

Let Tx, fix) ^ L'^( 7 x, Ex, fix) be the class of Levy processes having 

independent components. Let L'^’^(7x, Ex, fix) := L'^(7x, Ex, fix), d eN. 

For a Borel measure V on and ^ G [0, oo)'^, we dehne a Borel measure 
V 02 ; on where {Voz){A) := Ylt=i xiV{AnAd,i) for a Borel A C Here 
. 4 . 1,1 := 1^ and Ad,i ■= {x = {xi, ..., Xd)' G : Xm = 0 for m 7 ^ /}, for d >2, 
1 < I < d. Recalling (2.24), introduce := o and oi := o. When z G [0, 00), 


y e 


E G 


pdxd 


and V is a Borel measure on 




we set y Oi z := zy, 


T, Oi z := zE and V Oi z := zV. Recall (2.4)- (2.5). 


We collect some formulae for the associated canonical triplets of X T 
(see Theorem 30.1 in for the univariate subordination; see Theorem 3.3 
in [ 6 ] for the multivariate subordination). 


Lemma 4.2. Let fc G {1, d}. Let X ~ L‘^’^( 7 x, Ex, fix). Let T ~ H^) 

be independent of X. Then we have: 
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(a) X Of, T L'^(7xofcT, Hjcofer) with 

7xofcT = lx Ok Dt+ / X P(X(s) e dx) nT(ds), 

7[0,oo)J Vo<||a;||<l 

^XofcT = Ofc Dt , 

HxofcT (dx) = (IIx Ofc PT)(dx) + [ P{X{s) e dx) nr(ds). 

J [0,oo)J 

(b) For all t > 0 


P{{X Ofe T)(t) e dx} = f P{X{s) e dx) P{T{t) g ds). 

J [0,oo)*^ 


(c) If, in addition, Dt = 0 and /[oijfc cin'r(t) < oo then X T ~ 

Fy''(0,nxo,r). 


In Part (a) of Lemma 4.2 the dependence of T enters into the formulae 
in a linear fashion. As a result, if a process X is independently subordinated 
by a superposition of independent subordinators then it can be written (in 
distribution) as the sum of independent processes: 


Proposition 4.1. Let n > 1, k E {1, d} and X ~ ^x, hlx)- 

Let X,Ti,.. .,Tn be independent with Ti ~ SkiDxi, for l<l<n. Let 

T := Pk and Y := X T. Then we have: 

(a) T ~ Sk{DT, Ht) with Dt = Ei=i Dt, and Ht = Eh ■ 

(b) Y ~ L'^(7y, Sy, Ily) with 7y = Eh'yx^kTu = Ei=i^x°kTi and 
ny = Eh nxofcT, ■ 

(c) If Xi,... ,Xn are independent eopies of X, also being independent of 
Ti,...,T„, thenY = EhXiOkTi. 

(d) If, in addition, both X]r=i /[o pj '^nr;(^) < oo and EhDr, = 0, 

then Y ~ FV^i^O, Ily) and X T; ~ FV‘^{0, IIxo,,Ti) for all 1 < I < d. 


Proof, (a) is well known, but can alternatively be deduced from the Laplace 
transformation, (b) follows from Part (a), owing to Part (a) of Lemma 4.2 
(c) follows from Part (b). (d) follows from Part (a) as an implication of 
Part (c) of Lemma 4.2, □ 
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4.3 Proofs for Subsection 12.1 

Proof of Theorem 


2.1 


(a) Let F = 5oT ~ where 

T,B are independent with T ~ GGG^{a,T) and B ~ BM‘^{fr,T). Observe 


that (2.16) extends to A G C with 3fJA > 0. This follows from Schwarz’s 
principle of reflection: the proof of Theorem 24.11 of |S21 can be adapted to 


our situation. Let 6* G M'’* and set Xe := |||6*||| — i (/i, 0) such that = 


exp{—tXg). Now (2.20) follows from (2.16) via conditioning on T(f): 


E[exp{i{9,Yt))] = E[e = exp { - faAg - f / log[(x +A 6 ))/x] T(da:)} . 

J (0,oo) 


Here the right hand-side matches the formulae in (2.20) 


(b) (2.21) is shown in [21] (his Proposition 3.3), whereas (2.22) holds as with 

F.. TT7, ... .^d\ 

dy 

(see 1 ^, his Equation (4)) we have 


gY{s,r) = r'^ in (|2.22[) (r > 0, s G §'’*) and any Borel set A C 


nv(4) = 


dn 


Y 


dy 


{y) ^y = 


ls% Jo 


^ dr ds . 
dy 


□ 


Part (a) of Theorem 2.1 


Proof of Theorem 2.2. (a) We omit the proof as it is similar to the proof of 


(b) We decompose T into a superposition of independent subordinators T = 
^jc{i d} where T® := al, and 


T/:= lc,(AT,)AT,, t > 0 , 0 ^ J C { 1 ,..., d} , 


(4.1) 


0<s<£ 


with Gj as in (2.25). Here AT{t) = T{t)—T{ t—) for f > 0. Also, I : [0, cx)) —)■ 


[0, cxd)) denotes the identity. By Proposition 


4.1 


we have F = Ejc{i,....d} 
where (F*^) is a family of independent Levy processes with = B og [al) 
and F-^ = H ~ 0, Hj) with ~ A'^(0, H;^) for J A 0. For J A 0 

we have r(C'j) = 0 = 0 ^ F-^ = 0 Hj = 0. 


To see (2.28), suppose det S > 0 and J 7 ^ 0 with T{Gj) > 0. Note 
tj 

S''’*( 0 , H^f) with, using its polar representation. 


dH^ = lcjnsl(s) 1 ( 0 , 00) (r) k{s,r) a{ds)dr/r , 
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where k{s,r) is the quantity in (2.17)-(2.18). In view of Lemma 4.2 


nj(dx) = / F(/iOt + e da:)n^(dt), 

Jcj 

where Z is a d-dimensional standard normal vector. 

As both det S > 0 and n^(Cj) > 0, Ily must be absolutely continuous 
with respect to ij, admitting the following density, for y G Vj, 


dn, 


POO roo 


di 


iy) = 


exp{-rT-l\\y-ryos\\l^^} 


J Jcjnsl Jo Jo r( 27 rr)#F 2 Yl^^j 


dr /C(s, dr) Q!(ds). 


Here we set ||a:|| j,, := XljeJ ^ ^ ^ ^ Expanding a square 

in the exponent yields 


]^\\y-ryos\\%, = ^Ibllllj,, - 


E 77 

j'eJ 


I J.s t 


so as to evaluate the interior dr-integral using the identity (2.12) for the 
modihed Bessel function K of the second kind. For ?/ G we get 


dn, 

dl.i 




(4.2) 


teJ 


X 


’ CjOSi J (0,oo) . 




ieJ 


#J/4 


X 


/W({2r+i: Il!/ll4.) • 


j&J 




where the RHS of (4.2) matches (2.28). In (2.28), observe 

r^'^(dn^/d£j)(rs) = 5 (j(s,r), r > 0, s G §|; fl Vj , 


where the RHS matches (2.29). This completes the proof of Part (b). □ 


4.4 Proofs for Subsection 2.2 


Proof of Proposition 2.1. (a) Let 0 < g < 1. Pick e: > 0 such that, for all 
r > 0, 


< ST-^ I r^-i e-^ dr < 1 A r-F 


(4.3) 
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By (2.17)~(2.18), we get from Fubini’s theorem and a simple substitution 
that 


'0<||2||<1 


|^||'?nT(d^) = 


[ [ r^||s||^e ^'^/C(s,dr)—Q!(ds) 

10<r<l J(0,oc) ^ 


'[0,oo)^ Jo 


' [0,oo)^ 


^<?-ie-lhlk dr r(dx). 

^Il-Q f e~^ dr T{dx). (4.4) 


In view of (2.15) and (4.3)-(4.4), T(da;)/||a;||'? is hnite if and only if 

/o<||z||<i (dz) is, completing the proof of (a). 

(b) Let p,t > 0. Pick £ > 0 such that, for all r > 0, 


t^T-P 1 


0<r<l 


< er-P 


•■r ' e ' dr < lo<r<ir " + l^jie' 


(4.5) 


Using similar arguments as in the proof of Part (a), we get from (2.17)-(2.18), 
Fubini’s theorem and a simple substitution that 



z\\PUT{dz) 





dr T{dx). 


(4.6) 


In view of (2.15) and (|4(5)-(4.6), we see that ij| 2 ||>i ikP nr(dz) is hnite if 
and only if /o<|| 3 ,||<i T(da;)/||a:||^ is, completing the proof of (b). □ 

Proof of Proposition 2.2. Let k = d. Let B where S is a 

diagonal matrix. Let Z = {Zi ..., Zd)' G be a standard normal vector, 
that is a vector with independent standard normal components. For s G 
[0, ooY introduce B*{s) := (S o sY^'^Z. By self-similarity oi B — pi we can 


write B{s) = pos + B*{s) = pos-f (Sos)^/^Z for all (but hxed) s G [0, ooY- 
For X = {xk)i<k<d £ and M = {rnki)i<k,i<d ^ the maximum 

norms are denoted by ||(r||oo = max^ |xfc| and ||M||oo = maxfc,/ \mki\, respec¬ 
tively. Let II • Hop be the operator norm of || ■ ||. The equivalence of norms in 
hnite dimensions applies to as well as we hnd a common constant 

Coo G [1, cx)) such that || • ||oo < C'ooH • || as well as || - Hop < (Pool! • i|oo- Further, 
with II • He denoting the Euclidean norm on there exists Ce G (0, cx)) such 
that H ■ 11 — ^e\\ ■ He- 
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(a) Let 0<g<2, and introduce a function h : [0, oo)^ —)• [0,1] by 


/.(.):= E[||i?(.)||n(o,i](||S(.)||)], .e[0,oo)^ 


‘(2.30)^(2.31)’: We get from the self-similarity oi B — fil that 


his) < Bl||B(s)|n = E[||,^os + (Eos)i/2z||’], 56 ( 0 , 00 )“ 


and, thus, with Ci := (\\fxp + ^[11X1/2^119]), 

h{s) < 2’''‘(||.05||J,||Mr+||.05‘““||2,B|||E‘'“Z||«|) 

< 2'>''‘Ci (II . 05||?„ ||f,||“ + II . os'/^ll^ BlIlE'/^Zin) 
= 2«''‘C^ (||5||^||Mr+||5||^2ii'|||E'/2zr]) 

< C'l(lo<||s||<l||s||““^ + l||s||>l||s||“) , s 6 [ 0 ,oo)f. 


As h is globally bounded by 1 , there thus exists C 2 G ( 0 , 00) such that 
h{s) < (72(11511^/^ A 1) for s e [0, oo)f. (2.31) follows from this and the finite¬ 
ness of ^| 5 ||<i ||s||''^^ nT(ds). (Recall ny(d|/) = P{B{s) e dy)IlT{ds) by 
Lemma 4.2[) The proof is completed by Part (a) of Proposition 2.1 


‘(2.31)^(2.30)’: Assume det X 7 ^ 0. For s G 
maxfc XfcfcSfc > 0 and P{Z G M^) = P{T}PZ G 
such that, since ||Z|| I ~ Xl = r(d/2,1/2), for any p G 


[0, oo)f observe ||Xs||oo = 
Rd) = p{B*{s) eRi) = 1 


P{ exp{p||R*(s)|U} \\B*{s)\\'^ lo<||B*(.)||<i) = 0} = P{\\B*{s)\\ > 1} 

< PmilPiClWEsU-^} < 1 . (4.7) 


For A := —X and s G [0, cx))f we get from Girsanov’s theorem that 
B‘(s) -^lOs = B‘{s) + (Eo5)A ~ £(B-(5)|exp {(A,S*(o)> - /||A|||„}P) , 
and, by self-similarity of i? — pi, 


h{s) = E[||/iOS + R*(s)||no<||;,os+B*(^)ll<l] 

= exp{ - ||Ai||^,/2} E[exp{-(A,R*(s))} ||5*(s)||'' 1o<||b*(^)II<i] 

> ||s||''/2exp{-||s||oo||A|||/2} ho(s/||s||) 

> ||s|r/2exp{-Goo||A|||/2}ho(5/||s||), (4.8) 
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for s G [ 0 , oo)^ with ||s|| < 1. (Note ||s||oo < C'oo||s|| < C^o for ||s|| < 1.) 
Here ho : — )■ [0, cxd) is defined by 

h„(s) := £[exp{-||A|U l|B-(s)|U} ||B-(s)r 1o<i|b.(.)II<.] ■ 


As -P(n '^kk'^k 7 ^ 0 ) = 1 , s H-)■ l(o,i)(||-B*(s)||) is lower semicontinuous on 
almost snrely (indicator fnnctions of open sets are lower semicontinnons). As 
a result, ho is itself lower semicontinuous by Fatou’s lemma with a compact 
domain Si. In particular, ho attains its global minimum at some sq G Si. 


Note ho(so) > 0 by (4.7) 

C3||s||''/21o<|p||<i for s G 
( 0 ,oo). The proof is completed by an application of Part (a 


To summarise, we get from (4.8) that h(s) > 
[0,cx))^ with C 3 := ho(so) exp{-Coo|| A|||/2} G 

of Proposi¬ 


tion 2.1, using similar arguments as in the proof of ’(2.30 )^(2.31)’. 


(b) Let t,p > 0. Recall Y = B o^T for independent B and T. 


‘(2.32)^(2.33)’: Set C 4 := As in the proof of Part (a) note 


B[iiBwri < c 4 {iisinipr+iisr''"B[iis‘'"zin} (4.9) 

< C 4 {£|iiE‘''"zf] + iisir(iipir+£|iiE''"zri)}(4.10) 


for s G [0, cx))'^. If /i = 0 then it follows from (|4.9|) that £'[||y(t)p] = 
E[E[\\B{T{t))\\P\T{t)]] < C 4 E[\\E^/^Z\\p] 4 F[||T(t)f “and the LHS is hnite 
provided i7[||T(t)||fo^] is. Otherwise, if /i 7 ^ 0, we get from ( |4.10[ ) that 

B[||F(()|n < C 4 {B[||E‘''=Zr]+£)|||T(()|n(||pr+£)|||E'-'=Zf|)}, 


and the LHS is hnite provided i7[||T(t)||^] is. In view of Part (b) of Proposi- 

>([23I|’. 

[ 0 ,oo) by g{s) := 


tion 


2.1 this completes the proof of ‘(2.30) = 

Dehne g 


‘(2.33)^(2.32)’: Suppose /i = 0. Dehne : Si —)■ _ 
F'[||i?*(s)||Pl||B.(s)||>i]. Employing similar arguments as in the proof of Part (a) 
we hnd Sq G S!^ such that inf§d g = miugd g = g{sQ) > 0 and, thus. 


ziiiB(ii)r] = BiiiB-(s)ir] = iii'ii’’''‘BiiiB*(viisii)ri > 9(so)iisr'7 


for s G [0, cx))i This extends to i7[||R(s)p] > ( 7 (so)||s||^/^ for s G [0, cx))'’*, 
including the origin. In particular, this implies the inequality £'[||y(s)||^] > 
g{so)E[\\T{t)\\P^‘^] by conditioning on T(t). In view of Part (b) of Proposi¬ 
tion | 2 . 1 [ this completes the proof for p = 0 . 
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Assume fXk 7 ^ 0. By the equivalence of norms, we have || ■ || < C^\\ ■ || 
for some C5 G (0, cxd) such that, for s G [0, cxd)*^, 

||s|| < C5 ||s||oo < C'slll/zrlloo ||/iOs||oo < (^ooCsIll/zilloo llhosll, 

miUfc |/ifc| 

and, with Cq := (CooC' 5 ||l//i||oo)^ 2 ^^P, using the self-similarity of -B — /xJ, 


||s|r < {C^C4l/^I\\^rE[\\^IOs + B*{s)-B*{sW] 

< Cs£|||B(s)r| + C6B1 ||B-(s)||»| < CsB|||B(s)r| + 

where C 7 := C‘^CeE[\\'E^/‘^Z\\P]. Thus, we can hud tq G (0, cx)) such that 
CeE[\\B{s)\\P] > ||sp — > ||s||^/2 for s G [ 0 , 00 )“^ with ||s|| > tq. 

To summarise, we have ||sp < 2C'8-E[||-B(s)||^] -|- for s G [0, 00 )'^ and, 
thus, E\\T{t)\\P < 2C8-B[||F(s)p] T-rg, by conditioning, completing the proof 
of 

(We omit the proof for k = 1 being similar but simpler.) □ 


4.5 Proofs for Subsection 12.41 


Proof of Proposition 2.3. (a) Let Cx := {0} U ([0, cxD)f\C>A). It is straight¬ 
forwardly checked that Cx is closed under taking convex combinations. For 
X G Ca we have ||xp < (A,x) < (FHA||£;||x|| and, thus, ||x|| < C'HAHe by the 
Cauchy-Schwarz inequality (|| -He denotes the Euclidean norm and C G (0, cxd) 
is any constant with || ■ ||e < C\\ ■ ||.) Thus, Ca is a bounded subset of In 
particular, Ca is a compact, as it is also a closed subset of Mf. 

Continuity of Sx is obvious. For x G C>a we have Hxp — (A,x) > 0 and 
thus Sx{x) G [0, oo)f, as desired. 

(b) Let A G M'^. We get from Fubini’s theorem and (2.17) that 


eA.^>n-r(dx) = 


,,r({A,s}-r)(^ 


lhll>i 


'( 0 ,oo) J 

roo 


ho, 00 )^ 


exp < — r 


X — 


X 


/C(s, dr)a(ds) 

Ah W r(dx) 

r 




(4.11) 

Consequently, if T([0, cx))f\(!7A) > 0 then A ^ Vt- For the remaining part, 
assume T([0, oo)f\Ox) = 0, and choose £ > 0 such that for all r > 0 


log (r) < £ 


■— < log (r) - 1 - e" 
r 


(4.12) 
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Note that 


exp{((A,a:) — ||a;|p)/||x||}T(da:) < sup x 


r(dx). (4.13) 


'Ox 


'Ox 


In (4.13) the right hand-side is hnite in view of (2.15). The proof of Part (a) 
is completed by combining (4.11), (4.12) and (4.13). 

(c) Suppose A G Vt- (By Part (b) we have T([0, oo)f\Ox) = 0.) In view of 
Part (a), S\ ■. 0\ ^ [0, cx))^ is Borel measurable. In particular, the image 
measure, denoted by 7 a = (T|C>a) ° of the restiction of T to 0 \ under 
S\ is a well-dehned Borel measure on 

With S\ as in (2.39), note that there is a constant C G (1, oo) such that, 
for X e Ox with ||iSa(x)|| > 1, 


\x\ 


= 1 + 


(A,x) 


< 1 - 


(A,x) 


JiSa(x)|| ||a:||2 — (A,x) ||x||2 —(A,x) 
and, thus, by the transformation theorem. 


< 1 + 


(A,x) 


X 


<C. (4.14) 


/ (1+log i|x||) A — 7A(dx) = / (1+log ||5a(x)||)A 
'[0,oo)^ ll^ll 


r(dx), 


'Ox 


<C (1 + log ||x||) A 

Jox 11^1 


ll‘5A(a;)|| 

r(dx)+ [ log- ||5A(x)||r(dx). (4.15) 


'Ox 


(To show the inequality, split Ox into {x E Ox ■ ||5a(x)|| < 1} U {x G (Pa : 
||‘5A(a^)|| > 1} and recall C* > 1 in (4.14).) In view of (2.15) and (2.40), the 
RHS in(4.15) is hnite, completing the proof. 

(d) Let A G t > 0. li k = d, then we have T = B o^T for independent T 
and B ~ BM‘^{ij,, S). Conditioning on T{t) yields 

Eexp{X,Y{t)) = Eexp{A,R(T(t))) = Eexp{(/ioT(t),A) + ^||A|||^y( 4 )} 
= Eexp{qx,d,T{t)) . (4.16) 

Otherwise, if fc = 1 then E exp {\,Y(t)) = E exp{qx^iT(t)}. In either way, 
this completes the proof of Part (d). □ 


Proof of Theorem 2.3. Let k = d, t > 0, \ E T>y- Let q := qx^ E 
in (2.41). As A G Py, we must have q E Vt by Part (d) of Proposition 2.3 


as 
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Let a = 0. Adapting arguments from the proof of Theorem 25.17 of |S2], 
e.g., we get from (2.16) that, for with q — E [ 0 ,cxd)'^, 


Eexp{z,T{t)) = exp < —t 


[0,oo)d* 


log 


|x||^ — {z, x) 


X 


r(dx)}. (4.17) 


Let Oq as in ( |2.38| ), but with A replaced by q. Observe T([0, oo)f\Oq) = 0, the 
latter by Part (b) of Proposition 2.3, Note that 5q(x)/||5q(x)|p = x/(||x|p — 
(g, x)) for X G Oq. Set /ia = /i + SA. 

Since (4.16) extends as well, we get from (4.17) that, still with a = 0, 


= exp i — t 


= exp 


'o, 


'o. 


log 

log 


|x|p - (g,x) - i {fixOx,e) + || 


I Sotc 


|x|P — (g, x) 


\Sq{x)r-i{^,xoSq{x),e) + l\\er 


l‘5,(x)i|2 


r(dx) 
r{dx 


Next, apply the transformation theorem to T\Oq and Sq : Oq ^ [0, cx))f to 
see that the RHS of the last display matches (2.27), but with a, /r, T replaced 
by 0,g,\,Tq, respectively, where Tq is the well-defined Thorin measure in 
Part (c) of Proposition 2.3, but with A replaced by q. 


V 


According to (2.27), if a 7 ^ 0 it is possible to decompose Y = B + Yo into 
independent B, Yq where B ~ BM^{fj.oa, Soa) and Yq ~ fi, S, T). 

Using the independence, the proof is completed for k = dhy noting that 


Eexp{X + ie,B{t)) /Eexp{X,B{t)) = exp {ti (0 ,/ia o a)--||0|||^^} , 6 e 


The proof of the remaining case, where fc = 1, is similar, but simpler. This 
completes the proof of the theorem. □ 
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